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Abstract 


A  general  method  is  presented  for  finding  asymptotic  solutions 
of  problems  in  vave -propagation.   The  method  is  applicable  to  linear 
symmetric -hyperbolic  partial  differential  equations  and  to  the  integro- 
differential  equations  for  the  electromagnetic  field  in  a  dispersive 
medium.   These  equations  may  involve  a  large  parameter  X.   In  the 
electromagnetic  case  X  is  a  characteristic  frequency  of  the  medium. 
The  parameter  may  also  appear  in  initial  data  or  in  the  source  terms 
of  the  equations,   in  a  variety  of  different  ways.   This  gives  rise 
to  a  variety  of  different  types  of  asymptotic  solutions.   The  expansion 
procedure  is  a  "ray  method",  i.e.,  all  the  functions  that  appear  in 
the  expansion  satisfy  ordinary  differential  equations  along  certain 
space-time  curves  called  rays.   In  general,  these  rays  do  not  lie  on 
characteristic  surfaces,  but  may,  for  example,  fill  out  the  interior 
of  a  characteristic  hypercone.   They  are  associated  with  an  appropriately 
defined  "group  velocity".   In  subsequent  papers  the "ray" method  deve- 
loped here  will  be  applied  to  acoustic    wave  propagation,  Cerenkov 
radiation,  transition  radiation,  and  other  phenomena  of  wave -propagation. 
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1.   Introduction 

In  recent  years  asymptotic  methods  have  been  developed  for  the 
solution  of  certain  boundary -value  and  intial-value  problems  for  linear  partial 
differential  equations.   These  problems  involve  a  parameter,  and  the 
methods  provide  one  or  more  terms  of  the  asymptotic  expansion,  say  for 
large  values  of  the  parameter,  of  the  solution  of  the  problem.   They  are 
often  applicable  to  problems  for  which  no  exact  solution  method  is  known, 
and  even  for  problems  which  can  be  solved  exactly  it  frequently  happens 
that  only  the  asymptotic  expansion  of  the  solution  is  sufficiently  simple 
to  be  useful  in  practical  applications.   Furthermore,  it  is  invariably 
true  that  the  methods  which  yield  the  asymptotic  expansion  directly  are 
very  much  simpler  than  the  procedure  which  involves  first  finding  the 
exact  solution  and  then  its  asymptotic  expansion. 

An  important  class  of  asymptotic  methods  is  characterized  by 
the  fact  that  certain  curves,  often  called  "rays",  play  a  central  role 
in  the  theory.   The  rays  are  of  fundamental  importance  because  all  of 
the  functions  which  make  up  the  various  terms  of  the  expansion  can  be 
shown  to  satisfy  ordinary  differential  equations  along  these  curves.   Often 
these  equations  can  be  solved  to  yield  explicit  formulas  for  the  asymptotic 
solution  of  a  given  problem. 

The  "ray  method"  has  been  extensively  developed,  primarily  by 
J.  B.  Keller  and  his  co-workers  at  New  York  University,  for  certain  linear 
elliptic  partial  differential  equations  involving  a  parameter.  When 
applied  to  the  reduced  wave  equation  or  the  time-reduced  form  of  Maxwell's 
equations  the  method  yields  a  "geometrical  theory  of  diffraction"  [7] 
which  generalizes  the  classical  theory  of  geometrical  optics.   The  extensive 
literature  on  the  asymptotic  theory  of  the  reduced  wave  equation  and  Max- 
well' s  equations  has  recently  been  unified  and  summarized  in  [12]  . 

More  recently  ray  methods  have  been  applied  to  linear  hyperbolic 
partial  differential  equations  [^,10,11,13]  such  as  the  (time -dependent) 
wave -equation 

u,  ,  -  c^(X)  Au  =  0.  (1) 
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This  equation  involves  no  parameter,  but  a  parameter; jnay  appear  in  the 
data  of  a  problem  or  in  a  source  term  which  is  not  indicated  in  (l). 
Other  hyperbolic  equations  such  as 

u,  ,  -  c^(X)Au  +  X^^(X)u  =  0  (2) 

"C"D 

may  involve  a  large  parameter  X. 

From  the  standpoint  of  the  asymptotic  method,  (l)  and  (2)  are 
the  simplest  representatives  of  two  important  classes  of  linear  hyperbolic 
equations,  namely  the  "dispersive"  and  "non-dispersive"  hyperbolic  equations 
that  conserve  energy.   Ray  methods  for  these  two  classes  have  several 
important  differences,  the  most  striking  being  the  nature  of  the  rays 
themselves.   In  both  cases,  the  rays  are  space -time  curves,  but  for  non- 
dispersive  equations  they  lie  on  "characteristic  hyper surfaces"  and  are 
called  "bicharacteristics".   For  dispersive  equations  they  do  not  lie  on 
characteristic  hypersurfaces,  and  may,  for  example,  fill  out  the  interior 
of  a  characteristic  hypercone.   Applications  of  the  ray  method  to  various 
problems  for  (l)  and  (2)  have  recently  been  treated  in  LllJ  and  [lOj. 

The  purpose  of  the  present  paper  is  to  study  the  application  of 
the  asymptotic  theory  to  equations  more  general  then  (l)  and  (2),  and  to 
develop  a  method  of  solution  which  is  applicable  to  a  large  variety  of 
problems  in  wave  propagation.   In  order  to  avoid  undue  complications  and 
to  maintain  close  contact  with  physical  problems  we  have  not  attempted  to 
find  the  most  general  class  of  linear  hyperbolic  equations  to  which  our 
method  is  applicable.   Most  of  the  hyperbolic  systems  of  equations  which 
are  Important  in  mathematical  physics  fall  in  the  class  of  "symmetric- 
hyperbolic"  equations.  We  shall  show,  in  section  2.3,  that  our  method 
applies  to  "asymptotically  conservative  symmetric-hyperbolic  equations." 
But  oiir  main  emphasis  in  this  paper  is  not  on  these  equations. 

Perhaps  the  most  important  classical  problems  in  wave  propagation 
are  concerned  with  electromagnetic  waves.   Tiie  electric  and  magnetic  fields 
satisfy  Maxwell's  equations  which,  for  free  space,  are  symmetric -hyperbolic 
equations.   But  for  dispersive  media  the  equations  are  not  symraetric- 
hyperbolic,  in  fact  they  are  not  even  differential  equations.   The  time- 
dependent  form  of  Maxwell's  equations  for  dispersive  media  are  integro- 
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differential  equations.   They  are  discussed  in  section  2.1.   In  order  to 
encompass  this  important  class  of  problems  we  have  formulated  our  theory 
in  terms  of  a  "general  system  of  equations"  vhich  includes,  as  special 
cases,  the  integro-differential  equations  of  electromagnetic  theory  as 
well  as  the  symmetric -hyperbolic  differential  equations.   The  "general 
system"  is  described  in  section  2.2. 

As  pointed  out  in  section  2.1,  in  the  electromagnetic  case,  the 
large  expansion  parameter  }v.  is  a  characteristic  frequency  of  the  mediinn. 
The  true  meaning  of  our  asymptotic  expansion  is  better  understood,  however, 
by  introducing  an  equivalent  dimensionless  parameter  X  =  X  a/c  »  1. 
Here  "a"  is  a  characteristic  dimensicn  of  the  problem  and  c  is  the  speed  of 
light. 

Let  us  indicate  briefly  the  steps  involved  in  the  asymptotic 
method  which  is  described  in  chapter  3.   For  problems  which  can  be  solved 
exactly  (such  as  those  discussed  in  chapter  k)    examination  of  the  asymptotic 
expansion  of  the  solution  shows  that  it  consists  of  a  sum  of  terms,  each 
of  which  is  an  asymptotic  power  series  in  X    involving  a  "phase  function" 
and  an  infinite  sequence  of  "amplitude  functions".  For  complex  problems, 
we  therefore  ass\jme  that  the  solution  is  also  given  by  a  stun  of  such  series. 
By  inserting  such  a  series  into  the  general  system  of  equations,  we  find  first 
the  the  phase  function  satisfies  a  first  order  partial  differential  equation, 
the  "dispersion  equation",  which  can  be  solved  by  the  method  of  character- 
istics.  The  characteristic  curves  are  the  rays  which  we  have  mentioned. 
The  characteristic  (ordinary)  differential  equations  which  define  these 
curves  will  be  called  "ray  equations".   By  the  method  of  characteristics, 
the  phase  function  can  be  obtained  by  integrating  a  simple  ordinary  differ- 
ential equation  along  the  rays.  We  also  find  that  the  amplitude  fimctions 
satisfy  ordinary  differential  equations  along  the  rays  and  that,  in  the 
most  important  cases,  these  "transport  equations"  can  be  solved.   In  order 
to  find  the  rays,  and  the  phase  and  amplitude  functions,  initial  conditions 
for  all  these  ordinary  differential  equations  are  required.   In  some  cases 
the  required  intial  conditions  follow  directly  from  the  data  of  the  problem. 
In  others  the  initial  conditions  are  obtained  by  an  "indirect  method". 
Here  one  solves  a  "canonical  problem"  which  has  the  same  local  features 


as  the  given  problem.   It  is^  however,  sufficiently  simple  to  be  solved 
exactly.   The  required  intial  conditions  for  the  given  problem  are  ob- 
tained by  examination  of  the  asymptotic  expansion  of  the  solution  of  the 
canonical  problem. 

The  indirect  method  was  originated  by  J.  B,  Keller  for  the 
solution  of  "diffraction"  problems  for  the  reduced  wave  equation.   In  these 
problems  "diffracted  rays"  are  produced  when  ordinary  rays  are  incident 
upon  edges  or  corners  of  boundary  surfaces,  or  are  tangent  to  smooth 
surfaces.   Then  the  initial  conditions  for  the  diffracted  rays  are 
determined  from  a  canonical  problem  involving  a  boundary  with  local  features 
resembling  those  of  the  given  problem.   These  canonical  problems  are  often 
quite  difficult  because  of  the  presence  of  the  boundaries.   In  our  work 
we  consider  systems  of  equations  much  more  complex  than  the  reduced  wave 
equation.   However  we  are  concerned  here  with  problems  of  "radiation" 
from  sources,  which  are  much  simpler  than  the  diffraction  problems.*  Thus 
the  canonical  problems  we  must  consider  are  not  complicated  by  the  presence 
of  boundaries.   They  are  problems  for  the  general  system  of  equations  with 
constant  coefficients,  in  infinite  space;  whereas  the  original  problem 
may  involve  variable  coefficients  and  boundaries.   Chapter  h   is  devoted 
to  the  solution  of  the  canonical  problems  which  we  require. 

In  chapter  5  we  find  the  asymptotic  solution  of  two  representative 
radiation  problems.   These  two  problems  are  intended  only  to  illustrate 
the  asymptotic  method.   A  much  larger  variety  of  problems  for  (l)  and  (2) 
are  discussed  in  [HJ  and  [^10  J .   Other  problems  for  more  general  systems 
of  equations  are  under  investigation  by  the  asymptotic  method.   In  particular, 
these  involve  intial -boundary  value  problems  for  symmetric -hyperbolic 
equations  [_k    ']  ,  and  the  very  interesting  problems  of  "Cerenkov  radiation" 
[5  ~j    and  "transition  radiation"  [^   2~\    for  the  electromagnetic  field 
equations  and  other  equations. 


^Problems  involving  reflection  and  refraction  by  surfaces  and  interfaces 
are  considered  in  [4  ].   A  diffraction  problem  for  (2)  is  described  in 
\_1    J .   Diffraction  problems  for  more  general  equations  have  not  yet  been 
investigated. 


In  chapter  6  ve  investigate  the  special  features  of  our  asymptotic 
solutions  which  appear  when  the  general  systems  of  equations  is  non- 
dispersive.   In  particular  we  find  that  we  may  then  superpose  asymptotic 
solutions  to  obtain  one  which  involves  an  arbitrary  "wave-form"  function. 
Such  "progressing  wave"  solutions  have  been  studied  for  a  variety  of 
hyperbolic  equations  [3,13'!.   In  [^   llj  progressing  waves  are  applied  to 
problems  for  the  wave  equation  (l)  involving  amplitude -modulated  waves, 
radiation  from  moving  sources,  and  reflection  by  moving  surfaces.   Although 
specific  problems  are  not  treated  in  chapter  6,  the  discussion  there 
illuminates  the  relationship  between  the  dispersive  and  non-dispersive 
cases. 

The  introduction  of  the  general  system  of  equations  in  chapter 
2  was  largely  motivated  by  the  electromagnetic  field  equations  for  dis- 
persive media.   In  chapter  7  we  return  to  those  equations  in  order  to 
demonstrate  that  the  general  theory  which  we  have  developed  is  indeed 
applicable  to  the  electromagnetic  case.   In  section  7-2  the  solution  of 
the  transport  equations  for  inhomogeneous  isotropic  media  is  analyzed  in 
detail.   This  leads  to  a  discussion  of  the  rotation  of  the  electric  polariza- 
tion vector. 

Throughout  the  paper  certain  conditions  are  imposed  on  the 
equations  considered  in  order  that  the  asymptotic  method  may  be  applicable. 
These  conditions  are  listed  for  ready  reference  in  appendix  G.   A  brief 
summary  of  notational  conventions  is  given  in  appendix  F. 
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2.   The  General  System  of  Equations 

2.1  The  electromagnetic  field  equations  for  dispersive  media 
In  Gaussian  units.  Maxwell's  equations  take  the  fonn 

|-D-CVXH=  -knJ  I-  B  +  C  V  X  E  =  0,  (l) 

V  .  D  =  4itp^  ^  ■  1  =  0-  (2) 

The  source  terms  p(t,X)  and  j(t,X)  must  satisfy  the  continuity  equation 

1^  +v.j  =  o.  (3) 

From  (l)  and  (3)  we  see  that 

I  (V.B)  =  0,  I  (V.D  -  W)  =  0.  (k) 

Hence  the  equations  (2)  are  automatically  satisfied  if  they  are  satisfied 
at  any  given  time  t.   For  a  non-dispersive  medium  equations  (l)  and  (2) 
are  supplemented  by  the  constitutive  equations 

D  =  e  E,  B  =  11  H.  (5) 

Here  e,   M-  are  scalar  or  matrix  functions  of  X  =  (x,  ,  x  ,  x_)  which 
characterize  the  mediiJin.   It  is  well  known  that  the  solution  of  (1,2,5) 
with  given  source  terms  p,  J  satisfying  (3)  is  uniquely  determined  by 
the  initial  values  E(0,X),  H(0,X)  of  the  electric  and  magnetic  fields, 
say  at  t  =  0.   If  the  initial  conditions  satisfy  (2)  then  those  equations 
are  satisfied  for  all  t  >  0. 

The  fourier  (time)  transforms  of  E,  H,  D,  B,  p,  J  satisfy  the 
"reduced  equations"  which  are  given  by  (l),  (2),  (3)  and  (5)  with  -^  replaced 
by  -i(jJ.   Here  10  is  the  transformation  variable.   Dispersive  media  are  then 

A 

characterized  by  the  fact  that  e  and  |j.  are  functions  of  lj  as  well  as  X. 
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It  is  important  to  note  that  then  the  constitutive   equations  for  the  time- 
dependent  fields  become  convolution  integral  relations.   If,  for  con- 
venience,  we  replace  E,  H  by  U^  ,  U_,  D,   B  by  V-,  ,  V_,  and  e,  p.  by  e  ,  e„ 
then  for  dispersive  media  the  constituitive  equations  for  the  transformed 
fields  are  given  by 

V  (w,  X)  =  e  (^,X)U  (w,X),         V  =  1,2.  (6) 

It  then  follows  from  elementary  properties  of  fourier  transforms  that  the 
time -dependent  fields  are  related  by 

V  ft,X)  =  [    f  (t,X)U  (t-T,X)dT;     V  =  1,2.  (7) 

The  functions  f  ft.X)  are,  of  course,  fourier  transforms  of  the  dielectric 
permeability  e,  =  e  and  the  magnetic  permeability   e  =  n: 

a  (u,  X)  =  fe^""^   f  (t,X)dt,    V  =   1,2.  (8) 

These  equations,  which  are  thoroughly  discussed  in  [  8  J,  have  an  interesting 
physical  interpretation:  We  see  from  (7)  that  at  each  point  X,  the  fields 
V  depend  on  the  fields  U  only  at  earlier  times,  provided  we  assume  that 

the  real  functions  f  satisfy  the  causality  condition 

V  

f  (t,X)  =  0  f or   t  <  0,  (9) 

This  equation  has  important  consequences  for  the  analytical  properties  of 
the  functions  e  which  are  discussed  in  [81. 

In  Gaussian  units  e  and  e  are  dimensionless.   Since  w  is  a 
frequency,  e^  must  be  a  ftmction  of  w/x  where  X  is  a  characteristic 
frequency  of  the  medium.'^  Thus  we  may  set 

ei,(^,X)  =  %(^>^)i  ^  =   w/x.  (10) 


^Specific  examples  of  £,  are  discussed  in  appendix  D. 
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Then  we  may  re -write  (8)  in  the  form 


f  (t,X)  =  I-   /  e"^^'"*  ?;(co,X)dc.;         v  =  1,2.  (ll) 

The  time -dependent  form  of  the  electromagnetic  field  equations 
for  dispersive  media  are  now  given  by 

^-c  VxU2  =  -W,       ^  -fc  Vxy^  =  0.  (12) 


V  (t,X)  =  /  f  fT,X)  U  (t-T,X)dT;  V  =  1,2.       (13) 

Here  we  have  omitted  equations  (2)  since,  as  we  have  seen,  those  equations 
will  be  satisfied  automatically  if  they  are  satisfied  at  a  given  time.   It 
is  clear  that  the  initial  values  of  the  fields  E  =  U-,  and  H  =  U^  do  not 
uniquely  deteraiine  the  solution  of  this  system  of  equations,  since  the 
functions  V  depend  on  the  values  of  the  functions  U,  for  t  <  0.   In  the 
absence  of  any  existence  and  uniqueness  theorems  we  assume  that  an  appropriate 
problem  for  this  system  of  equations  is  as  follows:   Source  functions  p  and 
J  satisfying  (3)  are  given  for  all  t  and  X.   Prior  to  some  initial  time, 
say  t  =  0,  fields  V  and  U  satisfying  (2),  (12),  and  (13)  are  given.   We 
are  then  required  to  find  the  values  of  V  and  U  for  t  >  0.   In  the 
simplest  case,  the  sources  and  fields  are  identically  zero,  for  t  <  0, 
and  the  sources  are  "switched  on"  at  t  =  0  or  later. 

Much  of  our  work  in  the  later  sections  will  be  concerned  with 
the  solution  of  such  a  problem,  asymptotically  for  X  ->«.   (The  functions 
f  and  £  depend  on  the  parameter  X.     A  typical  example  is  discussed 


_9-  ^'-'^ 

in  appendix  D. )  Since  1  has  the  dimensions  of  a  frequency,  the 

physical  meaning  of  our  expansion  will  be  clarified  if  we  can  find  an 

equivalent  dimensionless  expansion  parameter,  X  .   In  order  to  do  this  we 

define  a  characteristic  length,  a.   We  then  introduce  dimensionless  space 

and  time  variables  X  =  X/a ,  t  =  ct/a,  and  set 

~o  ^^'     '      o  '     ' 

Ivo^^o'ifo^  =  Yv^^'^);  U,„(t^,X^)  =  y^(t,X),  f,„(t^,X^)  =  |  f,(t,X),  x^  =  c  x/a, 

^vo^-'V    =   ^v^'^'^)^      io^^o'^o)  =  -c  ^S^'^l-  (1^) 

Then  (ll),  (12)  and  (I3)  become 

X    P      -iX  wt 

f   (t  ,X  )  =  ;r^   /   e   °  °  e   (co,X  )dco;     v  =  1,2;  (I5) 

vo  o'^o    2n     /  vo   ''~o   '  '    '  \    -^  1 


^^10  ^^20 

xT^  -  V  X  u  ^  =  -Wj  ,     3-^  +  V  X  n  „  =  0;              (16) 

dt  o   '«20      '"0-'     dt      o   -10    '                ^   ' 

o  o 


Ivo^^o'^o^  =  j   ^vo^^o'L^Hvo^^o-^o'^o^^-^o'      ^  =  1'2;       (I7) 
where 

X  =  Xa/c.  (18) 

It  is  not  difficult  to  see  that  if  we  find  the  asymptotic  expansion  for 
X  ^  00  of  a  solution  of  (l6,17)  the  result  is  equivalent,  under  the  above 
transformation,  to  the  asymptotic  expansion  for  X  -^  00  of  the  corresponding 
solution  of  (12, 13).   We  conclude  that  the  correct  interpretation  of  our 


*For  problems  with  boundaries,  we  may  choose  a  typical  boundary  dimension 

as  the  characteristic  length,  a.   For  problems  involving  inhomogeneous 

e'(co,X) 
media  we  may  use  an  average  value  of  the  functions  I  .   Then  e 

Mv£;kx)i 

changes  by  a  sm.all  fraction  of  itself  over  distances  small  compared  to  a. 
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expansion  is  that  it  is  valid  for  1  «  X  =  Xa/c. 

Before  proceeding  to  oixr  discussion  of  the  asymptotic  solution 
of  (12,13),  ■we  shall  rewrite  that  system  of  equations  in  matrix  form. 
The  matrix  notation  has  several  advantages.   It  is  not  only  compact  and 
suggestive,  but  more  important,  it  will  enable  us  to  include  in  our  general 
discussion  a  large  set  of  equations,  of  which  (12,13)  is  a  special  case. 
In  particular,  we  will  show  that  this  set  includes  the  important  class  of 
(weakly  dissipative)  symmetric  hyperbolic  partial  differential  equations. 
The  notation  used  is  explained  in  appendix  F.   We  first  introduce  column 
vectors  u,  v,  f  with  6  components,  defined  as  follows: 

u  =  (U^,U  )  =  (E,H)  =  (E  ,E  ,E  H  ,H  ,H  ),   v  =  (V  ,V  ) ,   f  =  (-4nJ,0). 

Corresponding  to  any  3-vector  Z  =  (Zt,z^,z_),  we  introduce  a  3  >  3  matrix 

wv  J.         i£         3 

(Z)    given  by 

0  -z^  z^ 


(Z)    = 


-Zn 


^1 


(20) 


Then  if  V  is  any  3-vector,   in  column  form,    (Z)V  =  Z  X  V.      We  now  set 

12        3 
K  =  (k   ,k   ,k   )    and  define  three   6x6  matrices     A   ,  A   ,   A     by  the  equation 

w>         123 


v=l 


kA'  = 

V 


0       -c(K) 
c(K)     0 


(21) 


(19) 


Here  we  have  used  matrix  block  notation.   (Since  k, ,  k  ,  k  are  arbitrary 
we  may,  for  example,  choose  k  =  1  and  k  =  k  =  0.   Then  (20)  and  (21)  give 
A  ,  etc.)  We  now  note  that  the  matrices  A  are  constant  and 
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A% 


=     Z 


^ 


v=l 


1      ^ 

V  =  l  V 


aIu 


0  -c(V) 

c(V)  0 


Hi 


u. 


-cV  X  n 


cV  X  u. 


(22) 


hence   (12)   takes  the  form 


V,    +  / A     u 


v=l 


X 


(23) 


Next  we   set 


F(t,X) 


f^(t,X)  0 

0 


^2^^'l^ 


f  (oo,X)    = 


e'^(oj,X)  0 

0 


^2^'^'^^ 


{2k) 


For  anisotropic  media  e  and  f  are  3-^3  matrices.   In  the  isotropic  case 
they  are  scalars,  and  in  (2U)  they  multiply  the  3  ><  3  identity  matrix  I^ 
which  we  omit.   Now  (ll)  and  (13)  become 


and 


F(t,X)  =  ^  /  e-i^^*  ?(c.,X)da), 


v(t,X)  =  /  F(T,X)u(t  -  T,X)dT. 


If  we  increase  the  number  of  space  dimensions  from  3  to  n,  then  X 
and  (23)  becomes 


n 
+  11  A% 


v^  +  /    A  u   =  f . 

V=l         V 


(25) 


(26) 


\  ^2 }  •  • '  }  ^_, ) 


(27) 


This  is  the  matrix  form  of  (II-I3) 


(2.2) 
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2.2  The  general  system  of  equations 

In  the  succeeding  sections  of  this  paper  we  shall  discuss 

systems  of  equations  of  the  form  (2.1.26,  T)  where  u,  v  are  column  vectors 

V         fP  V 

of  dimension  m  and  A  ,  F  and  (2,  are  m  X  m  matrices.   The  matrices  A  are 

hermitian  and,  in  general,  may  be  functions  of  X_.   In  this  paper  we 
require  that  A  and  G  ^e  smooth  functions'*  of  X  in  the  infinite  X-space. 
If  appropriate  boundary  conditions  are  specified  on  boundary  surfaces, 
or  if  A  and  G    are  discontinuous  at  interior  surfaces  called  interfaces, 
our  methods  must  be  modified  by  introducing  reflected,  transmitted,  and 
diffracted  waves.   These  features  are  discussed  in  [l,  2,  k,  7,   11,  12j . 
In  addition  the  function  G=  (5  ((J,X;X)  is  required  to  satisfy 
five  "basic  conditions"  which  are  listed  in  appendix  G.   The  significance 
of  these  conditions  (and  the  other  "special  conditions")  are  best 
appreciated  by  examining  the  role  they  play  in  later  sections.   For  the 
electromagnetic  field  equations  in  isotropic  media,  examples  of  the 
dielectric  and  magnetic  permeability  functions  satisfying  these  con- 
ditions are  discussed  in  appendix  D. 


*A  "smooth"  function  posesses  sufficient  continuity  or  differentiability 
properties  to  justify  the  operations  performed  with  it. 
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2.3  S.ymmetrlc-h.yperbolic  partial  differential  equations 

In  this  section  we  shall  show  that,  for  an  appropriate  choice 
of  (^  ,  consistent  with  the  basic  conditions,  (2.1.26,  7)  reduce  to  a 
system  of  symmetric -hyperbolic  partial  differential  equations.  We  first 
note  that  (2.1.25,  6)  yield 

v(t,X)  =  /'F(t-s,X)u(s,X)ds  =  l^jTe-^^'^^^-^^cf  (co,X)u(s,X)dajds,    (l) 
and  the  fourier  integral  theorem  implies  that 

u(t,X)  =  ^  r^-±Mt-s)   ^(3  x)dc.ds.  (2) 

We  now  set 

where  A(X) ,  C  (X)  and  C  (X)  are  herraitian,  B(X)  is  anti-hermitian 
(i.e.  B  =  -B) ,  A  is  positive -definite,  and  C,  is  non -negative.*  It 
follows  then  from  (l)  and  (2)  that 

V  =  Au  +  XBu  +  Cu,  (h) 


hence  (2.1.27)  becomes 


n 


Au^  +  /    A^  u   +  XBu  +  Cu  =  f .  (5) 

^^X,  '"•X       '***    ■^'^        "^ 

V=l  V 


*The  matrix  Q  is  positive  definite  (or  non -negative)  if  (x,Qg)  >  0 


[or  (x,i^)  >  Oj  for  every  non-zero  vector  x- 


-1^- 


(2.3) 


If  C  =  0,  the  symmetric -hyperbolic**  system  (5)  conserves  energy  because 
B  is  anti-hermitian.   We  call  (5)  asymptotically  conservative  because  the 
large  parameter  X  does  not  multiply  the  dissipative  term  Cu.   The 
asymptotic  theory  of  such  systems  is  discussed  in  [^J . 

Since  A  =^(w^)=Ais  easy  to  see  that  (3)  satisfies  conditions 
1  2  and  h   of  appendix  G.   To  see  that  condition  5  is  satisfied  we  note 
that  G  =  '^+  (iX)~  C  -  toA,  where  ^=k^A  -  IB  +  x~  C  .  ^    is  hermitian 
because  B  is  anti-hermitian.   If  w  is  a  solution  of  det  G  =  0  then  there 
is  a  null  eigenvector  r  such  that  Gr  =  0.   The  equation  (r,  Gr;)=  0 
now   leads  to 

(r,(^r)   ^  (£,  C^r) 

and  condition  5  follows  from  the  fact  that  A  is  positive  definite  and  C^ 
is  non-negative.   Condition  3  imposes  an  additional  restriction  on  the 

V 

matrices  A  ,  A  and  B. 


**Actually  (5)  is  symmetric -hyperbolic  only  if  the  coefficient  matrices 
are  real  and  satisfy  the  conditions  we  have  imposed.   (Then  A  and  A  are 
symmetric  matrices  and  A  is  positive  definite.) 
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3.   As.ymptotic  Solution  of  the  General  System  of  Equations 

3.1  The  asymptotic  expansion 

In  this  section  we  shall  develop  a  method  for  obtaining  asymptotic 
solutions  of  the  general  system  of  equations  (2.1.26),  (2.1.27)  with  f  =  0. 
For  convenience  we  introduce  the  summation  convention  ( see  appendix  F) , 
and  frequently  we  shall  not  indicate  explicitly  the  dependence  of  various 
functions  on  the  variable  X.   Thus  we  may  write  the  general  system  of 
equations  in  the  form 

v_^  +  A^  u   =0,  (1) 

Wvt  l~X 

V 

v(t)  =  /  F(t)  u(t-T)dT.  (2) 

WV  J  -- 

From  (2.L25)  we  see  that 

F(t)  =  XI^(Xt),    ?f(y)  =  ^    f    e-^"^^  ?(c.)du).  (3) 

If  we  introduce  the  transformation  Xt  =  y  then  (2)  becomes 

v(t)  =  /  7^(y)u(t  -y/x)  dy.  (4) 

We  now  assume  that  (l)  and  (2)  have  solutions  given  by  asymptotic 

power  series  of  the  form 

00 

u(t,X)  =  e^^^^'l^  Yl     (i^)"""  \   (t.X).  (5) 

m=0 

We  shall  determine  the  phase  function  s  and  the  amplitude  function  z  =  z 

by  inserting  (5)  in  (l)  and  (2).*  We  introduce  the  convenient  notation** 

k^  =  s^  (t,X),      CO  =  -s^(t,X),  (6) 

V 


■'f-Presumably,  the  lower  order  terms,  Z-,,z  , ...,  can  be  obtained  by  a  similar 
procedure. 

**Thi3  notation  is  motivated  by  the  consideration  of  "plane-wave"  solutions. 

(See  the  footnote  to  section  3-2) 
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and  expand  u(t-y/x)  for  large  X.   In  order  to  do  so  we  first  note  that 


i\s(t-y/x)  _  iXs(t)+iwy 


2 
y  s 


2iA 


^*  +  0(X"^) 


]■ 


(7) 


and 


00 

m=0       "~  w~  L 


Hence 


■% 


u(t-yA)  =  e^^"^*)-'^"^  ]  z(t)  +  (iX)-^|z,(t)-iyz,(t)  -  |y^s^^z(t)|  +  oCx'^)' 


iv-^  --y-t' 


From  (3)  we  see  that 


(9) 


i(ca)  =  J^e^'^^;^(y)dy,   g^  =  jT  iye^'^^ff  (y)dy,   g^  =  -/y^e^^y^(y)dy,  (lO) 


and  we  apply  condition  1: 

!'((.)  =  l'(co)  -  (iX)-^^(c.O  +0(X-2) 
Then  (i|),  (9),  (10),  and  (11)  yield 


(11) 


v(t)  =  e^'^'^^^^Si.).   f  (iX)-l[g(.)z^-^(.)z  -g^z^  -f  I  s^^^^z]  -.  0(X-2)j  .   (12 


By  differentiating  (12)  and  (5)  we  obtain 

v,(t)  =  e    J  -iXcj  gz  +  ^z,  -wgz^  +  to^z  +  oj  g 

+  o(x-i)|  , 


cj-'t   2  tt^ww'-   ^  oj  t- 


(13) 


^] 


and 


.V      iXs 
A  u   =  e 


■  iXk^A^'z  +  Fa'^z^  +  k/^zJ  +  0(x"^)|  .        -{Ik) 


We  may  now  insert  (I3)  and  (1^4-)  in  (l)  and  equate  to  zero  the 
coefficients  of  X  and  X  .   The  resulting  equations  can  be  written  more 
simply  if  we  introduce  the  matrices 
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G  =  k^A^"  -  ox?  =  k^ix)    -  03  E(co,X)  (15) 

and 

A°  =   (oo^)^  =   -V  (16) 

Then 

A°  =  A°  co^   =   -s^^n^g        +  2g  1  ,  fl7") 

t  oj      t  tt  [_       coco  *- wj  '  V-^''' 

and  (1),  (13),  and  (ll|)  yield 

G  z  =  0  ..  (18) 


and 


G^i  +  a\  +  a\  +  I  A°  z  +  c^z  =  0.  (19) 

V 


From  (l8)  and  (19)  we  shall  derive  the  equations  which  determine  s  and  z. 
In  order  to  do  so  we  must  first  examine  some  properties  of  G. 
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3.2  The  eigenvectors  of  G 

From  (3.1.18)  we  see  that  the  vector  function  z  can  be  non- 
trivial  only  if 

det  G  =  detTk  A^(X)  -  wg(co,X)l  =  0.  (l) 

For  real  values  of  u,  X,  and  K,  this  equation  defines  a  functional  relation 
between  these  quantities  which  we  call  the  dispersion  relation.-^  The 
same  functional  relation  can  also  be  expressed  in  the  form 

CO  =  h(K,X),  (2) 

Then,  in  general,  h  is  a  multiple -valued  function  of  K  and  X,  defined  by  (l). 

For  each  value  to  =  h,  the  hermitian  matrix  G  is  singular, 
therefore  there  exists  a  positive  integer  q  <  m,    and  q  linearly  in- 
dependent null  eigenvectors  r  ,.  .  .  ,  r  of  G,  such  that 

Gr"^  =0;         j  =  l,...,q.  (3) 

The  integer  q  is  the  nullity  (or  the  multiplicity  of  the  zero  eigenvalue) 
of  G.   The  nullity  corresponding  to  each  value  of  (2)  may  be  different.   In 
our  work  we  shall  require  that,  for  each  value  of  (2),  q  is  independent  of 
K,  so  that  we  may  differentiate  (3)  with  respect  to  k^.   This  requirement 
is  a  consequence  of  condition  3  of   appendix  G. 

By  condition  2  the  matrix  A  is  positive  definite  for  all  real  w. 
Therefore  for  each  value  of  (2)  we  may  orthonormalize**  the  null  eigenvectors 
by  the  condition 

*In  case  A^  and  5  are  independent  of  X,  and  cf  =  ^,  the  dispersion  relation 
can  be  obtained  by  assuming  "plane  wave  solutions"  of  (3.1.1,2)  of  the  form 
u  =  z  exp[iX(k  X  -  wt)l  where  z,  k  ,  and  w  are  constants. 
**This  orthonomalization  may,  for  example,  be  accomplished  by  the  Gram- 
Schmidt  method  using  the  inner  product  |x,jj  =  {x,Aj^). 
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(r^,  aV)  =  5   ;     i,j  =1,...,  q.  {h) 

The  inner  product  used  here  is  defined  in  appendix  F.   5. .  is  the  Kroenecker 
symbol. 

If  we  differentiate  G  =  k  A  -  w  £  with  respect  to  k  and  recall 
that  A  =  x~  (w  £)  we  obtain 

\     =  a"   -   g^A°;      V  =  1,...,  n;  (5) 

V 


where 


=1  =I^S^5)'      v=l,...,n.  (6) 

V       V 


Hence  differentiation  of  (3)  yields 

G  r^  +  (a""  -  g  A°)r<^  =  0,'  (?) 

V 

and  from  {k) 

{t^,A^t^)   =  g     b      ;       i,j=l,...,q;  v=l,...,n.         (8) 

This  basic  identity  will  be  used  repeatedly  in  our  further  work.   From 

(k)   we  see  that  it  holds  also  for  v  =  0,  if  we  define 

So  =  1-  (9) 

The  eigenvectors  r  ,  •••,   r  are  not  completely  determined  by 

(3)  and  (if),  for  if  we  set 


/\±       \         V 
r   =  /    B.   r  : 

v=l 


i  =  l,..,q,  (10) 


/^l 


where  (3-  )  is  any  q-dimensional  unitary  matrix,  the  r   still  satisfy  (3)j 
(h) ,   and  even  (8).   Let  us  now  suppose  that  condition  6  of  appendix  G  is 
satisfied.   Then  if  we  set  a.  .  =  (r  ,L0"Or  )  the  q-diraensional  matrix  (a.  .) 
is  hermitian.   If  we  should  now  take  ((3.  )  to  be  the  unitary  matrix  which 
diagonalizes  (a.  ),  then  we  would  find  tha  the  new  eigenvectors  r  would 
satisfy  the  identity 

(r^,   ^^r^)    =r].?>^y  i,j=l,...,q.  (ll) 

Therefore,  if  condition  6  holds,  we  may  choose  eigenvectors  r  ,...,r^^  which 
satisfy  (3),  (h),    (8),  and  (ll). 
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3.3  The  dispersion  equation  and  the  ray  equations 

The  dispersion  relation,  introduced  in  the  last  section,  is  a 
functional  relation  between  w,  K,  and  X.   Since  w  =  -s_,  and  k  =  s 

V 

(3.2.1  or  2)  is  also  a  first  order  partial  differential  equation  for  the 

phase  function  s(t,X).   ¥e  call  that  equation  the  dispersion  equation. 

Some  readers  may  find  it  convenient  to  view  the  dispersion  equation  in  the 

form  (3.2.2)   as  a  Hamilton -Jacobi  equation.   In  any  case  it  may  be  solved 

by  the  "method  of  characteristics"  [  3  J •   Thus  we  introduce  the  characteristic 

(Hamilton's)  equations 

dx        >,      dk      ^, 

V        dh     V     dh  ,  /-.  N 

dt    ^v   dk  '   dt      dx  '         '        '  ^    ^ 

V  V 

An  immediate  consequence  of  these  equations,  and  the  fact  that  the  hamilton- 

ian  h  is  independent  of  t,  is 

du  _  dli  _  3h  V   5h  V  _  ,    s 

dt  ~  dt  ~  Sc^  dt    cK^  dt  ~  ^    ' 


From  (1)   we   see  that 

djc, 
dt   "  'Sc     d^  "*"  ^ 


ds        5s     V        Ss  _    /7  (^\ 


where  X  is  the  lagrangian 

V 

(1)  is  a  system  of  2n  first  order  ordinary  differential  equations  for 
the  2n  functions  x  (t),  k  (t).   Each  solution  of  this  system  defines 
a  CTorve  Ft,  X(t)l  =  Ft,  x^,...,x  1  in  space-time  which  we  shall  call 
a  ray.* 


■^Sometimes  the  term  "ray"  is  used  for  the  projections  X(t)  of  these  curves 
into  X-space. 
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If  initial  values  for  X  and  K  are  known  then, in  many  cases,  the 
ray  equations  (l)  may  be  solved  to  determine  a  ray:  Along  this  ray  we 
see  from  (2)  that  w  is  constant  and  from  (3)  that  s  may  be  obtained  by 
integration.   In  this  process  it  is  clear  that  we  require  initial  values 
for  X,  K,  oj,  and  s.   Let  us  examine  how  such  initial  values  are  obtained: 

¥e  shall  find  that  in  typical  problems  the  values  of  s  are  known 
on  some  manifold  "%  .   Thus,  for  example,  "{^  may  be  an  n-dimensional  hyper- 
surface  given  parametrically  in  terms  of  n  parameters    $  =  ( (f),  , . . .  ,  (|)  )  by 

X  =  X^($),     t  =  t($).  (5) 

Then  s  is  given  by 

s[t(£),  XJo)]  =  Sj$)  (6) 

where  s  is  a  known  function.   By  differentiating  (6)  we  see  that 
^        ^     Ss 

and  then  for  each  $   (7)  and  the  dispersion  relation  provide  n  +  1  equations 
for  the  determination  of  the  initial  values  w  ,  k   ,  ...,  k^^  of  the  n  +  1 
quantities  oj,  k,  , . . .  ,  k  . 

But,  in  general,  the  manifold -^Z may  be  of  dimension  r  <  n. 
Then,  for  the  p-urpose  of  our  asymptotic  method,  we  require  the  solution  s 
of  the  dispersion  equation  which  is  equal  to  s  on  "^ ,  and  for  which  the 
corresponding  rays  all  emanate  from  W with  increasing  t.   This  solution 
is,  in  general,  multiple -valued,  not  only  because  the  dispersion  relation 
defines  03  as  a  multiple -valued  function  of  K  and  X,  but  also  because, 
for  each  value  of  w,  more  than  one  ray  emanating  frcm'^  may  pass  through 
a  given  space-time  point  (t,X).   In  order  to  construct  this  solution  s 
we  note  that  (7)  and  the  dispersion  relation  with  X  =  X^(  $)    now  provide 


w>.      vvx,** 
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r  +  1  equations  relating  the  n  +  r  +  1  quantities 

cl)^,...,(i>,w,k^,...,k.  (8) 

T2_}        >  Tj,}      Q>     q2_>  'on 

We  may  then  introduce  any  n  of  these  quantities  or  any  n  functions 

7  .   =  7 -i^-,  >•  •  •  A    )'^   }^  -,  >•  • '  >^^     )  (9) 

of  these  quantities  as  independent  parameters,  provided  the  n  +  r  +  1 

equations  (9),  (7),  and  the  dispersion  relation,  can  be  solved  to  obtain 

the  quantities  (8)  as  functions  of  r  =  (7.^,.  . .  ,7^)  ■      The  functions  (9) 

can,  in  general,  be  so  chosen  that  (j)^, ,  . ,  ,cJ)^,co^,k^^, . . .  ,k^^  will  be 

single -valued  functions  of  r.   Then,  from  (5)  and  (6)  we  shall  have  the 

initial  values  X  ,  t,  s  ,  u  ,  K  given  as  single -valued  functions  of  r. 
v,.o'   '   o   o  ".O  ■~- 

Furthermore,  to  each  value  of  r,  there  corresponds  one  branch  of  the 
multiple -valued  function  h  for  which  w  (r)  =  h  K  (r),  X  (r)  .  With  this 
branch  of  h  we  may  solve  the  ray  equations  (l)  and  obtain  the  solution 

X  =  X(t;r),   K  =  K(t;r).  (lO) 

This  solution  satisfies  the  initial  conditions  X(T;r)  =  X  (r),  K(T;r)  = 


=  Kjr). 


The  parameters  r  lie  in  a  certain  parameter  space  (^  which  is 


determined  as  follows:   In  (5),  the  point  (t,X)  =  (t,X^)  ranges  over  the 
manifold'??  as   0=  (<)),.., ,(j)  )  ranges  over  some  space  ef  .   Tlien  f 
consists  of  all  values  r  =  (7,,...,  7  )  assumed  by  the  functions  (9)  when 
$   ranges  over  «f  and  co  ,  K  take  on  all  real  values  that  satisfy  (7) 

/**  o  '^  o 

and  the  dispersion  relation. 


(3.3) 
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For  r  In  P ,   and  t  >  T(r),  the  equation  (t,X)  =  |t,X(t;r)1 
defines  an  n-parameter  family  of  rays.   On  each  ray,  w  has  the  constant 
value 

cj  =  CO  (r)  =  hrK(t;r),x(t;r)].  (11) 

Wow  if  we  set 

g^(t;r)  =  h   [K(t;r),  X(t;r)1;     V  =  l,...,n  (12) 

"      V  -' 

and 

J?(t;r)  =  k  (t;r)g  (t;r)  -w(r);  (13) 

then  (3)  yields,  by  integration 

t 
s(t;r)  =  srt,x(t;r)l  -  s  (r)  +  f    i(t';r)dt'.  (li^) 

T(r) 

This  equation  provides  the  required  multiple -valued  solution  of  the 
dispersion  equation  which  is  equal  to  s  on  '^K  ,   and  for  which  the 
corresponding  rays  all  emanate  from  *%  with  increasing  t. 
Note:   In  this  and  succeeding  sections  the  symbol  f(t;r)  denotes  the 
value  of  the  function  f(t,X)  on  the  ray  rt,X(t;r)~I  at  time  t.  i.e. 

f(t;r)  =  f  rt,x(t;r)1  (15) 

For  fixed  t,   ^(^5^)  may  be  multiple -valued  although  f(t;r)  is  single- 
valued  because  more  than  one  ray  may  pass  through  the  point  (t,X). 


(3.^) 
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3.  h-     The  ray  transfomation 

For  each  fixed  value  of  t,  the  function  X  =  X(t;r)  defines  a 
transfonnation  from  the  parameter  space  fi   to  X-space.   The  jacohian 
j(t;r)  of  this  ray  transfonnation  is  defined  by 


j(t;r)  =  det 


S  7 


(1) 


The  jacobian  may  vanish  at  certain  space-time  points  called  caustic 
points.   The  locus  of  all  such  points  is  called  the  caustic  of  the  ray 
family. 

A  complete  study  of  j(t;r)  would  require  the  solution  of  the 
ray  equations  to  obtain  the  function  X(t:r).   However,  we  can  obtain 
valuable  information  about  the  behavior   of  j  near  ^  by  means  of  an 
expansion  of  X(t;r)  for  small  rt-T(r)l.   From  the  ray  equations  (3.3.I) 
we  obtain 


Then 


Six    Sx 

V  _   ov 

S7  ~  S7 


T   h^   +  O(t-T) 


(3) 


From  (1)  it  follows  that 

j(t)  =  j[T(r);r]  =  det 


3k 


ov 


W 


If  j(t)  =  0  for  some  value  of  T   then  the  point  (t,X)  =  ItjXCt;?) 
is  a  caustic  point,  and  we  find  from  (3)  that  there  is  a  positive  integer 
V  such  that 
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(3-M 


where 


j(t;r)  =  (t-T)' B'Ct;!;)  +  o[(t-T)^-'i] , 


TCTjr)    =     lim     (t-T)"""  jFt;?] 


(5) 


(6) 


The  quantity  j  is  non-zero  and  finite.   It  can  be  obtained  from  further 

terms  of  the  expansion  (3).   It  is  of  practical  importance  to  note  that 

V  and  j  can  be  obtained  from  the  ray  equations  (3.3-1)  and  the  initial 

conditions  X  (p),  K  ( P)  without  solving  the  ray  equations  explicitly- 
«-  o  ^   "^O  -~ 

In  general J  the  integer  v,   which  is  the  order  of  the  zero  of  j(t)  at 
t  =  i{v) ,   is  a  function  of  P.   If  j(t)  f  0   then  the  point  [T,X(T;r)]  is 
not  a  caustic  point.   However  (5)  and  (6)  remain  trivially  vaHd  with  v  =  0 
and  j(T;r)  =  j(t;P).   Thus  v  and  j  are  defined  at  all  points  on  every 
ray,   and  in  general  v  is  a  non -negative  integer. 

It  is  Interesting  to  note  that, since  g  =  h   , 


j(t)   =  det 


^1 


g. 


ol 


n 
ax 


on 


^1 


^7, 


ax 


ol 


ax 


on 


^ 


^. 


(7) 


This  equation  follows  from  (k)   by  simple  determinant  operations.   From  the 
ray  equations,  we  see  that  the  first  row-vector  of  (7)  is  tangent  to  the 
ray  [tjXftiPH  at  t  =  t(p).   The  remaining  row-vectors  are  tangent  vectors 
of  "??? .   If,  at  a  given  point,  the  n  tangent  vectors  are  linearly  dependent, 
then  either  that  point  is  a  singular  point  of  *%  or  7n   is  of  dimension  less 
than  n.  In  either  case  j(t)  vanishes  there,  hence  the  point  is  a  caustic 
point.   If  at  a  given  point,  the  n  tangent  vectors  are  linearly  independent, 
then  we  see  from  (T)  that  j(t)  vanishes  if  and  only  if  the  ray  at  that  point 


is  tangent  to  rM 
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3-5  The  transport  equations 

In  the  preceeding  section  we  have  seen  that  the  phase  function 
s  can  be  obtained  by  solving  a  system  of  ordinary  differential  equations , 
the  ray  equations.   In  this  section  we  shall  show  that  the  amplitude 
function  z  can  also  be  obtained  by  solving  a  system  of  ordinary  differential 
equations. 

From  (3-1-18)  we  see  that  the  vector  z  lies  in  the  null  space 

of  the  matrix  G,  hence  is  a  linear  combination  of  the  null  eigenvectors 

r  ,...,r  of  G.   Thus  there  exist  coefficients  a.,,..., a     such  that 
^   '        '-n.  1'    '  q 


(2) 


.=  tia     r"'.  (1) 

^   ;^  m^ 

f 

If  we  take  the  inner  product  of  (3- 1-19)  with  r  and  note  that 

t  t 

(r  ,  Gz^ )  =  (Gr  ,  z, )  =  (0,  z, )  =0,  we  obtain 

{r\   A^z   )  +  |(r\  A°  z)  +  (r^  c.  Jt)  z)  =  0;  (Yl 
J  \ J=0, 

Here,  for  convenience,  we  have  set 

x^  =  t.  (3) 

We  now  insert  (l)  in  (2)  and  simplify  the  resulting  equation  by  using 
the  basic  identity  (3.2.8).   Thus  we  obtain 

n   q 


J  J  ^J=0  m=l 

n  n         ,      J     - 

But  from  (3. 3.1)  ve  see  that  2_Sj(Vx  =^''lV/—  (^ X       ^=     ^  '     ' 

j=0      J      j=l      j 

hence  (k)   may  be  written  in  the  form 


(^) 
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.  , ^,   a  =  0;      I  =   l,,..,q  (5) 

m=l 


where 


''^   J=0 


=  Z_  (r  .  A^  r   )  +  -(r  ,  A  r  )  +  (r\  ojcOr  ); 

^,m  =  1,. . . ,q.       (6) 

(5)  is  a  system  of  q  first  order  ordinary  differential  equations  for  the 

coefficients  a-,,...,CT  which  determine  z.   These  equations  are  called 
1      q  fc** 

the  transport  equations.   The  solution  of  (5)  yields  the  values  of  cr,,...,CT  , 
hence  the  value  of  z,  along  a  ray.   In  the  next  four  sections  we  shall 
obtain  the  solution  of  the  transport  equations  in  several  important 
special  cases.   From  (3.2. 11)  we  see  that  if  condition  6  holds  the  last 
term  in  (6)  is 
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3.6  Homogeneous  media 

We  shall  say  that  our  general  system  of  equations  (2.1.26,7) 
describes  the  propagation  of  waves  in  a  homogeneous  medium  if  the  matrices 
A^,  g,  and^are  independent  of  X.   In  this  case  it  is  possible  to  show 
that  if  condition  6  holds  the  system  of  q  equations  (3' 5- 5)  uncouples.   In 


(1) 


(2) 


fact  it  follows  from  the  theorem  of  appendix  A  that 

n 

T,   =  ^  V.  G  +  T]  I5 ,;     V-  G  =  /~  (g  )   ; 
•—  _l  v=l      V 

and  we  obtain  the  uncoupled  equations 

dT -^  |(gv)x^  -^^^]^^  =  °'   t  =  i,...,q. 

For  homogeneous  media  h  and  g  are  independent  of  X.   It  then 
follows  from  the  ray  equations  that  the  functions  k  are  constant  on  each 
ray,  i.e. 

K'^'^'^t)  =  ^o^£^-  (3) 

Then,  of  course, 

w(t;r)  =  u)  (r)  =  hfK  (r)l,   gJt;r)  =  g  (r)  =  h  [k  (r)]  (k) 

and  the  solution  of  the  ray  equations  (33)  is  given  by 

\  =  \(^'E)  =  \o^E)  ■"  b-^^vj^vo^}:^'  (5) 

Thus  the  rays  form  an  n-parameter  family  of  straight  lines  in  space -time. 
From  (3- 3'^)  we  see  that  the  lagrangian  has  the  constant  value 
i(t;r)  =i  (r)  =  k   g   -  CO  (6)' 

^   '^       O^^  VO       VO      o  ^  ' 

on  each  ray,  hence  (3.3-1^)  becomes 

s(t;r)  =  s  fr)  +  rt-T(r)1^  fr).  (7) 

The  solution  of  (2)  is  given  in  (12)  of  appendix  C.   There  we  see 
that 
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^.(*'?  =  ^(^o'l:^ 


J(t^;r) 


-11/2 


j(t;r) 


exp    -(t-t^),^(r) 


t=i,, 


(8) 


Here   j   is  the  jacobian  defined  by  (S.l+.l),   and  t[     is  given  by 
n/r)    =   P(r),   ^^jSC'^Jr^r)];        r^r)    =r^a3     K  ); 


w     =  0)  ( r) ;         K    =  K  ( r) . 


(9) 


(8)  provides  the  value  of  a  at  any  point  on  the  ray  ['tjX(t;r)J  in  terms 
of  its  value  at  a  given  point  t  =  t  .   From  (3. 5.1)  and  (3.2.8)  we  see 
that  CT_  =  (r™,  A°z),   therefore  (3, 5-1), and  (8)  yield 

m 


m 


Z(t)   :. 


t 

m=l 


exp  ]  -(t-t  )ti  (a  (t  )r 
^  I     "■        o  'm  (  m^  o  - 


m= 


_exp{-.(t-t^),j[AA°z(tj]A    (10) 


Here,  for  the  sake  of  brevity,  we  do  not  indicate  the  dependence  on  the 
parameters  r.   If  condition  7  of  appendix  G  is  satisfied  the  functions 


■m 


m 


l,...,q;  are  all  equal.   Then  (10)  becomes 


z(t) 


j(t) 


1/2 


exp  I  -(t-t^)n|  z(t^), 


(11) 
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It  is  interesting  to  note  that  for  homogeneous  media 


-T7 =  (%o  -  ''  ^vo\     ■"  ^^vo\     ^'  (12) 

hence  the  jacobian  (3.^.l)  is  a  polynomial  in  t  with  coefficients  which 

are  functions  of  P.   The  degree  of  the  polynomial  is  at  most  n.   If 

T  (r),...,  T  (r)  are  the  roots  of  this  polynomial  then  the  factor 
1  i^       n  ■^ 

jC't  )/j('t)  which  appears  in  (8),  (10),  and  (11)  is  given  by 

j(t  )    (t  -T_,)...(t  -T  ] 
'^^    o    ^  o  1    ^  o  n 

r^)  •••  (t-T^) 
(If  the  degree  of  the  polynomial  is  less  than  n  then  some  of  the  roots 


iitr  =  (t-Tj  ...  (t-T )  •  (13) 


T.  in  (13)  are  infinite.)   The  points  t  =  T.(r)  are,  of  course,  caustic 
points.   At  those  points  (13)  becomes  infinite  and  our  formulas  for  the 
amplitude  z  are  invalid.   In  the  neighborhood  of  a  caustic  point,  a 
different  asymptotic  representation  of  the  solution  is  required.*  If  t 
and  t  are  separated  by  one  or  more  caustic  points  then  (13)  may  be  negative. 
In  such  cases  the  factor  j(t  )/j(t)l  '      in  (8),  (10),  and  (11)  is  ambiguous. 
The  ambiguity  may  be  removed  by  the  "phase-shift  rule"  which  will  be 
discussed  in  section  3.10. 


*See  [6] 
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3. 7  The  energy  density 

We  tiirn  now  to  a  consideration  of  the  transport  equations 
(3.5-5)  without  assuming  that  A  ,  §  and  (tQ  are  independent  of  X.   We 
first  note-  that  the  basic  identity  (3.2.8)  implies 

^^   '^  £t)  -"  ^^t'^   £  )  +  (£  '\l   )  =  (£  ^A  r  )^  =  i^J^   =  0;    (1) 
hence,  from  (3-5.6)  and  (3-2.8) 

*-*.-fl  [(eSaV   )    Mr'    ,  A^-)]    .    [r^^^.^*)r"■J  (2) 


tm         int       ^  L'-         -Xj,  --Xj^ 


=  Z:  [(rj,  a'  f^K    -  (L.  <   £")]  *  fc'.  "(»)£"] 


V=l  V  V 


n 


tm         -  -         ^^3_ 


where 


n 


V-G  = 


Z     (gJ.   .  (^) 


■'v'x 

V=l  V 


Next  we   introduce  the  energy  density  ftmction. 


=     (z,A°z)    =  o^^a^  f  2_  1  (5) 

which  is  discussed  in  appendix  B.   (The  second  equation  in  (5)  follows 
from  (3.5.1)  and  (3.2.8).)   The  transport  equations(3.  5-5)  now  yield 

g  =  -  a*  (t,  +  t  ,*)a   .  (6) 

dt      t   to    mt   m 

Comparison  of  (6)  and  (3)  shows  that  we  may  obtain  a  single  ordinary 

differential  equation  for  w  provided  conditions  J   and  8  are  satisfied 

(We  note  that  oondition  8    is  trivially  satisfied  if  the  matrices  A 

are  independent  of  X.   Furthermore,  conditions  7  and  8  are  obviously 

satisfied  if  the  nullity  is  q  =  1.) 
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Under  these   conditions,    (6)   now  becomes 

^     +   rv-G  +2(Re  Ti    +  p)]w  =  0  (7) 

cit  *-     ^■*» 

and  the  solution  of  this  equation  is  given  in  appendix  C.   Prom  (C.IO)  ve  see 

that 

J(t  ;r)      (         ^  ^ 

v(t;r)  =  w(t^;r)  -tjIt^  exp  j  -2  J[ReTi(t' ;£)  +  p(t' ;r)|df  -  .     (8) 

As  in  the  preceding  section  we  note  that  j  vanishes  at  caustic  points 
and  therefore  w  becomes  infinite  at  such  points.   Except  for  the  ex- 
ponential factor,  (8)  would  imply  that  w(t;r)j(t;r)  is  constant;  i.e. 
independent  of  t.   Let  dX(t)  be  an  element  of  space,  each  point  of  which 
moves  according  to  the  solution  X  =  X(t;r)  of  the  ray  equations.   Tlien 
dX(t)  =  j(t;r)dr,  and  the  element  of  energy  w(t;r)dX(t)  =  wCt;?) j (t;r)dr 
is  independent  of  t.   Thus  energy  is  conserved  in  dX(t).   Since  dX(t)  vanishes 
at  a  caustic  point  it  is  not  surprising  that  the  energy   density  becomes 
infinite  there. 

The  energy  density  w  is  by  definition  a  non-negative  real  quantity. 
However  j(t  )/j(t)  may  be  negative  if  t  and  t  are  separated  by  a  caustic 
point.   In  this  case  the  derivation  of  (8)  in  appendix  C  is  not  valid,   A 
generalization  of  (8),  valid  for  all  t  except  at  caustic  points  will  be 
given  in  section  3-10' 

In  this  section  we  have  obtained  a  formula  for  the  energy 
density  w,  but  not  for  the  amplitude  function  z.   For  many  purposes  this 
formula  suffices.   (We  recall,  e.g.,  that  often  optical  intensities, 
i.e.  average  energy  densities,  are  measiirable  whereas  the  corresponding 
fields  are  not.)  As  we  shall  see  in  the  next  two  sections,  further 
progress  can  be  made  in  the  important  special  cases  q  =  1  and  q  =  2. 
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3.8  Nullity  one 


If  q  =  1  then  conditions  7  and  8  are  trivially  satisfied  with 

n 


Ti(w)  =  []^,u}M  r]        and    3  =  "  o  ^  ^^' 


a"  r) 
^  v=l      ^v 


From  (3- 5. 5 J 6)  we  see  that  the  transport  equation  becomes 


g-fT   a  =  0 
dt  o 


(1) 


(2) 


where 


Li 

=  y__  (r,  A"^   r      )    +  -(r,   A°  r)    +  t] 


3=0 

Let  9=arg  cr.      Then 

I     I    i© 
a  =     a  e 


and   (2)   yields 

[dlcTl 


.  I   I  de 


+    \a\i      =   0. 
'     '    o 


(3) 


(h) 


(5) 


By  taking  the  imaginary  part  of  (5)  we  find  that 


dO 

dt  =  -^-^ 


,     hence  0(t)  =  Q(t^)  -  /   Im  T^(t')dt' 


(6) 


From  (3.7.5)  we  see  that  \a\    =  Jv ,     hence  (3-7.8)  yields 

1/2         r    t 


a(t)|  =  |a(t^) 


rJ(to)i 


J(t) 


A 


exp  /  -J  [Ee   T(t')+(3(t')]df> 


(7) 


Thus  from  (k) ,    (6),  and  (7) 


a(t)  =  a(t^) 


rJ(t  ) 


1/2 


j(t) 


f 


exp  /  -/  [Re  Ti(t')+P(t')-H  Im  T^(t')]dt'^ 
L  t 


(8) 


(3.8) 


Since  (3.1.17)  shows  that  A,  is  hermitian,  we  see  that  (r.  A,  r)  is  real. 
Hence,  from  (3) , 

Re  T]  +  i  Im  T  =  T)  +  1^  (9) 

where 


n 


Ylir, 


^  =  Im^  (r,  A^  r   ).  (lO) 

(if  A  ,    •••,   A  and  G  are  real  then  the  eigenvector  r  is  real  and  ^  =  0.) 
From  (3.5.1)  and  (3.2.8)  we  see  that 


z  =  or      and     a  =  (r,  A  z).  (11) 

Therefore  (8),  (9),  and  (ll)  yield 


z(t)  = 


j(t) 


,v°(  , 

o     o 

t 


exp  ^  -/  [n(t')+3(t')+iaf)>t'/  [r(t  )A°(t  )£(t  )]r(t).  (12) 


o 


Here,  for  the  sake  of  brevity  we  have  not  indicated  the  dependence  on  r  of 

1/2 

all  the  ftmctions  in  (12).  Again  the  ambiguity  of  the  factor  [j(t  )/j(t)]  ' 

may  be  removed  by  the  "phase-shift  rule"  of  section  (3. 10). 

We  note  that  the  eigenvector  r  is  determined  only  up  to  a  unitary 
factor  e   ,  where  a  is  an  arbitrary  real  number.   It  can  be  shown*,  however, 
that  (12)  is  independent  of  a. 


•^Suppose  r  =  e  j;.   Then  tJ  =  t),  3  =  p,  and       r    =iaer+er 

Hence  (r,  A'^r   )  =  (r,A'^r   )  +  ia  g,,  where  from  (3-2.8)  g.  =  (r,A'^r). 

It  follows  that  ?  =  ^  +  If  and   /  't(t')dt'  =  /  ^(f )dt' -ta(t)-a(t^) .   But 

t  t 

expr.i[,(t)-a(t^)])  [r(t^),A°(t^)z(t^)]r(t)  =  [^(t  J  ,  A°(t^)z(t^)]r(t) . 


I 
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(3.9) 


3-9  Nullity  two 

In  important  applications  of  our  theory,  such  as  the  electro- 
magnetic field  equations  in  isotropic  media,  the  nullity  q  of  the  matrix  G 
is  two.   In  such  cases  it  is  possible  to  solve  the  coupled  transport 

equations,  provided  we  Impose  conditions  6,  J ,   8,  and  9-   Then  we  may  choose 

1      2 
the  eigenvectors  r  and  r^  to  be  real  and  it  follows  from  (3. 5-6, 7)  that 

T   is  real.   In  this  case  (3-7.2)  implies  that 

^11  =  '^22  =  '^1  (^^ 

where 

T^  =  I  V.G  H-T  +  p.  (2) 

Similarly,  we  see  from  (3.7.2  )   that  t.„  +  Tp-,  =  0,  hence  (3.5.6  ,7)  yield 


'^12   "'^21   '^2 


(3) 


where 


Tp  =1_     (r\AJr2   )  +i(rl,A°r2)_  ^^^ 

The  transport  equations   (3.5«5)  now  take  the  form 

da  da 

IT  -^  -^1^1  ^  -^2^2  =  °^  dT  -  ^2^1  ^  ^1^2  =  ^-  ^5) 

These  equations  can  be  simplified  by  introducing  the  quantities 

P  =  +  /w        p^  =  a^/p,        pp  =  a^/p.  (6) 

Then,  since  p~  =  "w,    we  see  from  (3- 7-7)  that 

2p  1^  +  [V.G  +2(11+  p)]p^  =  0,  (7) 

hence  (2)  yields 

f  .  .,P  =  0.  (8) 

By  inserting  (6)  in  (5)  we  obtain 


-36- 
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dt" 


•^  V2 


0, 


dp. 
dt" 


T2P1  =  0. 


and  it  is  easily  seen  that  the  solution  of  (9)  is  given  by 
Pl(t)  =  p^(t^)  cos  &(t)  -  p^Ct^)  sin  &(t)  ] 


PjCt)  =  p3_(t^)  sin  5(t)  +  p^Ct^)  cos  5(t; 


(9) 


(10) 


where 


&(t)  =   /   T2(t')dt'. 


(11) 


(12) 


Here  p,  (t)  =  P-,(t;r)  is  the  value  of  p  at  the  point  t,X(t;r)1,  etc 
Wow,  from  (6)  and  (3' 7 .5)  we  see  that 

|3i|2  +  \^f     -     (v)-'  [|aj2^  la^l^j  =1, 

while  (3.5.1)  and  (6)  yield 

z(t)  =  a^(t)rl(t)  '+  a2(t)r2(t)  =  ^/^^{t)    [p^(t)rl(t)+P2(t)r2(tr|  .   (I3) 

Furthermore  (6),  (3.5.I),  (3-2.8)  and  (3.7-5)  imply  that 


/•  m  .0  X    /  m  ,0  ^ 
(r  ,A  z)    (r  ,A  z) 


Pm  = 


y(z,A°z) 


m  =  1,2. 


(1^) 


By  inserting  (3.7.8)  in  (I3)  and  by  using  (3.7.5)  we  obtain 

■(t  )1"^/^        ^ 

^—^  expj  -Jh(f)+p(t')]dt'j  [z(t^),A°(t^)z(t^)]^/2 


z(t) 


We 


X  [p^(t)rl(t)+P2(t)r2(t)J.    (15) 
now  set  p;(t)  =  [i,(t^),A°(t  )z(t  )]V2       ^  ^^^^^ 


o''^^  o--'     'm-  ■    '  ^  o'  ""m 
Then  p^  and  p^  satisfy  (10),  while  (ik)    and  (15)  yield 

m 
■\^'  A  ^^o^S^^o^J'       "^  =  1'2  (16) 

1/2      r  _. 


P;^V  =  &:"(V'A°^Vs(t  )]; 


and 


2(t) 


j(V 


exp(  -/ 


LTl(f  )+p(fj|dfy    [p^(t)rl(t)+p'(t)/(t)].        (17) 


^^s  if  2(^0^    ^^  Siven,    (16)    and   (10)   yield  p^(t)    and  p^(t),   and  then   (17) 
determines   z(t). 


(3.10) 
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3.10   Generalized  formulas  Tor  the  amplitude,  ^ 

In  sections  3. 6,  3.8,  and  3-9  we  have  derived  formulas  for  the 

amplitude  function  z(t;r)  in  terms  of  its  values  z(t  ;?)  at  some  point 

t  =  t  (r)  on  each  ray.   We  have  seen  that  z  becomes  infinite  at  caustic 
o  ~  •» 

points  because  the  jacobian  vanishes  there.   Thus  our  fonnulas  are  valid 
only  if  t  is  not  a  caustic  point. 

In  section  3.3  "we  derived  the  formula  (3- 3-1^)  foi"  "the  phase 
function  s  in  terms  of  its  values  on  the  manifold  /5^  ,  which  consists  of 
those  points  on  the  family  of  rays. for  which  t  =  T(r).   In  general,  these 
points  may  be  caustic  points. 

For  applications  of  our  theory  we  require  formulas  for  z  and  s 
in  terms  of  their  values  on  the  same  manifold  7^.   If  t  =  T(r)  is  not  a 
caustic  point,  we  can  merely  set  t  =  t  in  our  formulas  for  z.   However,  if 
t  =  T  is  a  caustic  point  z(T,r)  is  infinite  and  we  must  modify  our  formulas 
for  z  appropriately.   That  is  the  object  of  this  section. 

We  shall  make  use  of  the  analysis  of  the  behaviour  of  the  jacobian 
near  7fl  ,   discussed  in  section  3.^.   In  particular  we  shall  use  (3.^.6). 
Thus  if  we  let  t  -t  tend  to  zero  through  positive  values  in  (3.7.8)  we  find 
that 

w(t;r)  =  w(T;r)  jj^j^j  exp  j  -2j^eTi(t'  ;r)+  3(t'  ;r)]dt'  (  ,        (1) 

where 

w(T;r)  =   lim     (t-T)^w(t;r).  (2) 

(t-T)->0+ 

(The  limit  is  taken  through  positive  values  of  t  -t  because  if  t  <  t  <  t 
^  .  -^  o  o 

then  t  and  t  are  separated  by  a  caustic  point  and  (3-7.8)  is  not  valid). 
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(3.10) 


The  positive  f\inction  w  describes  the  asymptotic   behavior  of  w  as  t  -^  x. 
If  V  =  0  then  t  =  t  is  not  a  caustic  point,  w(T,r)  =  w(T;r),  j  =  J j  snd   (l) 
reduces  to  (3.7.8).   Equation  (l)  is  valid  under  conditions  7  and  8. 

If  we  let  (t  -t)  tend  to  zero  through  positive  values  in 
(3.9.10,11,16,17)  we  obtain 


z(t) 


j(t) 


11/2 


exp 


■/[Tl(f  )+p(f)Jdt'V  [^^(t)r^(t)-fP2(t)r2(t)];    (3) 


where 


P^(t)  =    lim     (t^-T)"/2p;(t^)  =  [V"(t),A°(t)z(t)];    m  =1,2;   (h) 

(t  -T)-^0-f 
^  O 


m  o 


z(t)  =  z(T;r)  =  liiTi    (t-T)""/^  z(t;r); 


(5) 


and 


3l(t)  =  P^(t)  cos  &(t)  -  ^^(t)  sin  5(t) 


};      &(t)  =  /  T2(t')dt'.  (6) 


p2(*)  =  Bi(^)  ^in  5(t)  -  ^^(t)  cos  &(t)  J  T. 

Equations  (3) >    i^)  >   and  (6)  are  valid  under  conditions  6,  7j  8,  and  9j  for  the 

case  of  nullity  q  =  2.   If  z(T;r)  is  known  then  (k)    gives  the  values  of 

P  (t),  (6)  gives  the  values  of  p  (t),  and  (3)  determines  z(t;r).   If  v  =  0 
m  m  -v^   wi 

then  J.  =  ^,  j  =  J}   and  the  formulas  reduce  to  those  of  section  3-9^  from 
which  they  were  derived. 
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(3.10) 


If  we  let  (t  -t)  tend  to  zero  through  positive  values  in  (3. 8.12) 


we  obtain 


z(t) 


HA 
j(t) 


1/2 


r  t 


exp  • 


[Ti(f)+3(t')^H^(f)]df>  [^(t),A  (T)z(T)]r(t),   (7) 


^   T 


where  'z   is  again  given  by  (5).   Equation  (7)  is  valid 

for  the  case  of  nullity  q  =  1,   It  provides  the  value  of  z(t)  directly 

in  terms  of  z(t),  and  reduces  to  (3.8.12)  if  v  =  0. 

For  the  case  of  homogneous  media,  we  may  let  ("t  -t)  tend  to 
zero  through  positive  values  in  (3. 6. 10,  11).   Then  if  condition  6  holds, 
we  obtain 


z(t) 


j(t) 


1I/2 


exp  \-{t-x)^ 


m=l 


m 


m 


.A°i(T)]: 


m 


(8) 


from  (3.6.10).   Here  t\    ,   r^ ,   and  A°  are  constant  on  a  ray  and  11^  is  given 
by  (3.2.11).   If  condition  7  is  also  satisfied  (3. 6. 11)  yields 


z(t)  = 


ilil 


-1/2 


jTt) 


exp  I  -(t-T)ri|  z(t) 


(9) 


In  (8)  and  (9)  z  is  again  given  by  (5).   If  v  =  0  (8)  and  (9)  reduce  to 
(3.6.10,  11). 

From  (3.J+.6)  we  see  that  j(t)  and  j(t)  have  the  same  sign 
provided  t  <  t  <  t' ,  where  t  =  T'(r)  is  the  location  of  the  next  caustic 
point  on  the  ray  X(t;r).   For  t  outside  this  interval  j(T)/j(t)  may 
be  negative.   We  assume  then  that  (l)  remains  valid  with  j/j  replaced 
'by  Ij/oI-  We  also  assume  that  the  other  formulas  of  this  section  may 
be  interpreted  according  to  the  following 


-^0-  (3.10) 


T3hase-shift*  rule: 

=  e'P 

[:i(t)1 

1/2 

^(t) 

[j(t)^ 

J(t) 

1/2 


For  T<t<T',P=0- 


p  is  constant  "between  caustic  points  and  decreases  by  —  when  a  ray 
passes  through  a  caustic  point  with  increasing  t.  Here  v  is  the  order 
of  the  zero  of  j(t)  at  that  caustic  point. 

The  phase -shift  rule  has  been  obtained  by  examining  the 
asymptotic  expansion  of  the  exact  solution  of  n-umerous  problems    \_6,   10  J 
It  forms  part  of  our  general  prescription  for  the  construction  of 

asymptotic  solutions.   (An  interesting  illustration  appears  at  the  end  of 
section  5-2) . 

It  is  interesting  to  note  that  for  homogeneous  media  w.K.r  .ri 

o 
and  A  are  constant  on  each  ray,  while  (3  =  0.   Then  it  is  easy  to 

demonstrate  that  (3)  reduces  to  (9)-   (The  demonstration  follows  from 

(3-6.1)  which  implies  that  t_  =  0,  and  from  (6)  which  implies  that 

P  and  p  are  constant).   It  is  also  easy  to  show  that  (7)  reduces  to 

(9)  if  condition  6  is  satisfied,  for  in  this  case  the  theorem  of  appendix 

A  show  that  t  is  real  and  (3.8.9)  implies  that  ^  =  0. 


*rhe  term  "phase-shift"  is  a  little  unfortunate  because  it  applies  to 
the  amplitude  function  z,   not  to  the  phase  function  s.   However  the  term 
is  in  common  use. 


.]+l.  (3-11) 


3.11  Asymptotic  Solutions 

Equations  (3,  h,   6),    (7),  (8),  and  (9)  of  the  preceding  section 
provide  formulas,  valid  under  different  conditions,  for  the  amplitude 
function  z(t;r)  in  terms  of  z(T;r).   In  addition  (3.3.1^1)  provides 
a  formula  for  s(t;r)  in  terms  of  the  value  s  (r)  of  s  on  ^ .  By  inserting 
these  results  into  (3.1.5),  we  obtain 

^^  e^^^^^'D  ^(t;r)  (1) 

for  the  leading  term  of  the  asymptotic  expansion  of  u.   This  equation, 
and  the  equation 

X  =  X(t;r),  (2) 

^lAl        kVkv         MW 

which  is  obtained  by  solving  the  ray  equations,  provide  a  parametric 
representation  of  u(t,X)  vith  n  parameters  T  =  (7, ,...,7  ).   For  a 
given  value  of(t,x)(2)  may  have  zero,  one  or  more  solutions  T.   It  is 
understood  that  for  each  (t,X)  (l)  is  to  be  summed  over  all  values  of  T 
which  satisfy  (2).   (In  geometric  terms,  we  sum  over  all  rays  of  the 
n-parameter  family  which  pass  through  a  given  space -time  point.) 

We  have  yet  to  answer  a  very  important  question:  How  do  we 
determine  the  functions  s  ( r)  and  z(T;r)  which  we  need  to  determine  s 
and  z?  It  is  possible  to  give  at  least  a  partial  answer  here. 

Obviously  s  and  z  must  be  determined  from  the  data  of  a 
suitable  problem  for  our  system  of  equations  (2.1.26,  27).   Often  these 
data  consist  of  the  values  of  the  source  function  f .   In  problems  involving 
boundaries,  the  boundary  values  also  provide  data.   In  initial  value  problems, 
which  are  appropriate  for  systems  such  as  those  described  in  section  2.3, 
the  initial  conditions  also  are  part  of  the  data  of  the  problem. 


(3.11) 
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In  some  cases,  the  values  of  s  and  z  can  be  determined  directly 

O       '^ 

from  the  data.   This  is  true,  e.g.  for  initial-value  problems  with 
"oscillatory  initial  data"  discussed  in  \^1,    k,   9,    lOj . 

In  other  cases,  such  as  those  which  we  shall  discuss  in  chapter 
^^  "M   is  of  dimension  r  <  n  and  consists  of  caustic  points.   The  values  of 
s  on  ■%  can  be  deduced  from  the  data  of  the  problem,  but  z  is  infinite 
on'h   and  the  values  of  z  must  be  obtained  by  an  indirect  method.   In  this 
method,  we  complete  the  construction  of  the  solution  (1,2)  with  the  function 
z(T;r)  unspecified.   Then  we  specialize  our  problem  to  one  with  constant 
coefficients,  i.e.  with  A  and  g  (to)  independent  of  X,  in  unbounded  X-space. 
(The  constant  values  chosen  for  A  and  E    {^)   are  the  values  of  A  (X)  and 
5'(w,X)  at  the  point  X  =  X(T,r).)  Thus  for  each  different  value  of  X(T;r) 
we  obtain  a  problem  with  constant  coefficients.   This  canonical  problem'^ 
is,  however,  much  simpler  than  the  original  problem;  and  it  can  be  solved 
exactly.   (This  is  done  in  chapter  h.)      Then  the  asymptotic  expansion  of 
the  solution  of  the  canonical  problem  is  found  to  be  identical  to  the  result 
obtained  above  (and  specialized  to  constant  coefficients),  except  that  'z(T;r) 
is  now  given  explicitly.   If  we  assume  that  z_  depends  only  on  the  local 
value  of  the  coefficients,  the  z   obtained  for  the  canonical  problem  can 
be  used  to  complete  the  construction  of  the  asymptotic  solution  of  the 
original  problem.   The  details  of  the  indirect  method  are  best  understood 
by  examining  it  in  specific  applications  such  as  those  of  chapter  5  and 
[l,  2,  h,    5,  7,  10,  11,  12 J. 


\  *rhe  term  "canonical  problem",  and  the  idea  of  the  "indirect  method"  were 

Introduced  by  J.  B.  Keller  in  studying  asymptotic  solutions  of  the  reduced 
wave  equation  and  the , phenomenon  of  diffraction  [6,  7,    12]. 


h.      Exact  solution  of  the  system  of  equations  with 
constant  coefficients 

U.l  The  fourier  integral  representation 

In  chapter  3  we  presented  a  method  for  obtaining  asymptotic 

solutions  of  the  system  of  equations  (2.1.26,  7)-   With  constant  coefficients, 

that  system  becomes 
n 
V  +2_  A%^  =  f(t,X),  (1) 


v(t,X)  =  /  F(T)u(t-T,X)  dr.  (2) 

Here  u  and  v  are  column  vectors  of  dimension  m,  X  =  (x,,...,  x  ),  and 


F(t)  =1^  j     e-i^"*  gH    dc..  (3) 

The  constant   mxm  matrices  A  are  hermitian,  and  the  mXm  matrix  6  {^) 

satisfies  conditions  I-5.   As  explained  in  section  2.1,  we  consider  the 

following  problem  for  (1,2): 

The  source  function  f is  identicially  zero  for  t  <  0.   It  is  required 

to  find  a  function  u(t,X)  which  satisfies  (1,2)  and  the  "initial  conditions", 
u(t,X)  =0        for  t  <  0.  (k) 

In  order  to  solve  (l,2,U)  we  set  K  =  (k^,...k  )  and  introduce  the 

fourier  transform 

u(oj,K)  =  /  exp  j  -iA(k  x  -  wt))  u(t,X)dt  dX.  (5) 


(U.l) 


Then 

u(t,X)  =  (^)     /   exp  fiX(k^x;-wt)]u(w,K)dLodK.  (6) 

Similarly  we  introduce  the  foixrier  transforms  v(uj,K)  and  f(w,K).   5('^) 
is,  of  course,  the  one -dimensional  fourier  transform  of  F(t).   Because  of 
the  convolution  integral  relation  (2),  it  is  easy  to  show  that 

v((J,K)  =  J(oj)u(u,K).  (7) 

From  (1),  we  find  that 

-iXcjv  +  iXk  A  u  =  f ,  (o; 

hence,  from  (7), 

(k^A"  -  wi)u  =  (iX)"^f.  (9) 

We  now  introduce  the  matrices 

G  =  k  A^  -  w  €(w),  G  =  k  A^-ojg(aj)  =  G+(iX)""^aj<2J+0(x'^) .  (lO) 
(Here  we  have  used  condition  1.)   Then  (9)  yields 

Gu  =  (iX)"^f;      u  =  (iX)-^  G"^f.  (11) 

Thus,  from  (6)  we  obtain  the  representation 

u(t,X)  =  (iX)"-^  (1^^     fexp   {iX(k^x^-wt)|  G"^((J,K)f(cj,K)da)dK.    (12) 

We  assume  that  f(t,X)  is  sufficiently  well-behaved  as  |xl  -»oo  to  justify 
the  interchanges  of  integrations  we  have  performed.  We  also  impose  the 
following  temporary  condition*  on  f : 

f  =  0     for  t  >  t^  >  0.  (13) 

This  condition,  and  the  fact  that  f  =  0  for  t  <  0,  imply  that 
(for  each  K)  f  is  analytic  in  the  entire  co-plane.  On  the  other  hand 
conditions  h   and  5  imply  that  G   is  analytic  in  the  half -plane  Im  w  >  0. 


*The  condition  is  removed  in  section  h.'y. 


45-  ^ 


In  (12)  the  integrations  with  respect  to  k^ , — ,k  are  taken  over  the 
real  axis.   We  choose  the  contour  in  the  co-plane  to  lie  just  above  the 
real  axis.   Then  for  t  <  0  ve  can  close  the  contour  (by  an  infinite 
semi-circle)  in  the  upper  half -plane,  and  since  the  integrand  has  no  sin- 
gularities there  the   resulting  integral  is  zero.   T^us  (k)   is  satisfied. 
Since,  by  construction,  (l)  and  (2)  are  satisfied,  (12)  is  the  solution 
of  the  problem  (1,2,4). 

k.2     Residue  evaluation  of  the  integral  representation.  ■ 

For  t  >  0  (and  for  each  value  of  K)  we  may  close  the  contour 
in  the  lower  half  of  the  u-plane  and  evaluate  the  to  integral  in  (U.1.12) 
as  a  sum  of  residues  at  the  roots  co  =  h(K)  of  det  G(w,K)  =  0.  Fran 
(U.1.10)  we  see  that  for  X  ->  oo  the  roots  Ti  tend  to  the  roots  w  =  h  of 
det  G(w,K}  =  0.   If  any  roots  u)  =  h  are  not  real  they  must,  by  condition 
5,  lie  in  the  half -plane  Im  w  <  0.   It  is  easy 'to  see  that  the  residue 
contributions  frcm  the  corresponding  roots  h  will  be  exponentially  small 
for  X  ^00,  and  will  not  contribute  to  the  leading  term  of  the  asymptotic 
expansion  of  u.  Therefore  we  need  find  the  residue  contributions  only  from 
those  roots  h  which  are  perturbations  of  real  roots  h.  '- 

We  assume  that  ina neighborhood  of  each  such  root,  G(u))  =  G(to,K) 

has  a  complete  set  of  eigenvectors  r  (co)  such  that 

G(a))rJ(a))  =  7J(w)rJ(u));    j=l,...,m.  (l) 

Therefore  we  may  introduce   the  expansion 

m 


f(u)) 


=   Yl  b  re)  ?J(u,),  (2) 


J=i     ' 


and 

m 


S-^((.)f(a))    =  Yl   %A^)/P(^)]lH^)  (3) 

.1=1    •-  "^  -• 
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At  least  one  eigenvalue  7"^(w)  of  G  vanishes  at  each  root  co  =  h  of  det  G  =  0. 
Suppose  7"^(w)  vanishes  at  the  point  w  =  h..   In  section  U.5  we  shall  find 
that  7'^(oj)  has  a  simple  zero  at  h.,  i.e. 

7J(h.)  f  0*  W 

.V 

It  follows  that  the  residue  contribution  to  (4.1.12)  from  the  point  u  =  h . 

J 

(and  from  the  j   eigenvalue)   is  given  by 

(-2rti)(iX)-^(V2n)''"'^  ^exp|iX(k^x^-h^(K)t)j  [b_.  (h_.)/f  ^  (hj)]rJ(h  .  )dK.   (5) 
Hence , 

Jti(*>$)  -^(iJ)  V  ^  ^^  (i^(  Vv-h.t^^  [-^j  (h.  )/;j(h.  )]rJ(h.  )dK.      (6) 

Here  we  sum  over  every  zero  h.  (of  every  eigenvalue  7  (w)  of  G)  which  is 
a  perturbation  of  a  real  zero  h  of  an  eigenvalue  of  G. 

In  order  to  simplify  (6)  further  we  must  examine  the  eigenvalue 
problem,  first  for  G  and  then  for  G.   This  will  be  done  in  the  next  two 
sections. 


*This  follows  from  (4.5-3)  and  {h.^.h)   for  i  =  j. 


■hi-  (^.3) 


4.3  The  eigenvalue  problem  for  G 


For  real  w  and  IC  =  (k-,,...jk  ),   let  7(w,K)  =  "/((ji)   be  an  eigenvalue 


of  multiplicity  q  of  G  =  k  A  -ajG(Lo).   Then  there  exist  q  linearly  independent 


V 


-oj  G  (lo) 
eigenvectors  r^((j^,K)    such  that 

Gr-^  =  7r^       j  =  l,...,q.  (l) 

By  differentiating  (l)  with  respect  to  w  we  find  that* 

-A°r"^  +  Gr"^  =  7  r"^  +  jr^ .  (2) 

Then,  since  (r^,  Gr"^)  =  (Gr^,  r"^)  =  /(r^.r"^)  we  see  from  (2)  that 

(r\  aV)  =  -7      (/,  r^);    i,j  =1,...,  q.  (3) 

By  condition  2,   we  may  orthonormalize  the  eigenvectors  with  respect  to  A  ; 
thus 

-7  (r^^r-^)  =  (r^,  aV)  =  6..;        i,j=l,...,q.       (k) 

If  condition  6  holds,  then  as  in  section  3-2  we  can  choose  the  eigenvectors 
so  that  they  satisfy  not  only  (l)  and  (k)   but  also 

[r^,  u)c2T((^)^J]  =  T^(c^)6.j;       i,j=l,...,q.  (5) 

Of  course,  if  co  =  h(K)  is  a  zero  of  7(oj)  then  at  cj  =  h 

GrJ  =0,  j  =  l,...,q.  (6) 


*Here  we  use  condition  3- 
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k.h     Perturbation  of  the  eigenvalue  problem 

Because  of  (Ji.l.lO)  we  assume  that  the  perturbed  matrix  G 
has  eigenvalues  of  the  form  .  _  :   ._:-.. 

7^(0))  =  r(w)  +  0(X"^),         j=l,...,ci,  (1) 

and  that  the  zeros  w  =  h .  of  ^"^(u)  are  perturbations  of  the  zeros  w  =  h  of  7 ; 

J 

h.  =  h  +  (iX)"-HiJ  +  0{\~^)  ,  j    =   l,...,q.  (2) 

J 

We  also  assume  that  the  null  eigenvectors  of  G(h.)  are  perturbations  of  the 

J 

null  eigenvectors  of  G(h) ,  i.e. 

rJ(h  )  =  rJ(h)  +  (iX)-^rj  +  0{X-^) ,        '  (3) 


where 


J 
and  {k)   yields 


G(h)r'^(h)  =  0     and     G(h .  )r,'^(h . )  =  0.  (k) 

J  J 

Now,  since  G  =  -A°'  and  G  =  G  +  (iX)"  co^(w)  +  0(x"  ), 

G(h.)  =  G(h)  -  (iX)-\J  A°(h)  +  (iX)-\^(h)  +  0(X-2),      (5) 


[h^(h)  -  hj  A°(h)]rJ(h)  +  G(h)r^  =  0.    ,  .,  (6) 


Thus 


[^^(h),  hg(h)rJ(h)]  =  hJfr^Ch),  A°(h)rJ(h)],  (?) 

and  we  see  from  (^•3«^^5)  that* 


h^  =  n  (h);      3   =   l,...,q.  (8) 


J 


*In  the  usual  perturbation  solution  of  an  eigenvalue  problem,  the  leading  term 
of  the  perturbed  eigenvector  ^  is,  in  general,  a  linear  combination  of  the 
unperturbed  eigenvectors  r^  ,  i  =  l,...,q.   We  have  succeeded  in  solving  the 
pertiirbed  problem  with  the  simpler  assumption  (3)  only  because  the  eigenvectors 
r  have  already  been  chosen  in  such  a  way  that  they  satisfy  (^•3«5). 


I4..5  As.ymptotlc  expansion  of  the  exact  solution 

We  are  now  in  a  position  to  obtain  the  leading  term  of  the 

asymptotic  expansion  of  (i|.2.6)  for  X  ^00,   We  first  use  (1|.2.2)  and 

(li..i4..3)    to  obtain 

m 
f(h.)  ~  zib.  (h.)  rJ(h).  (1) 

Then  (U.3.i|-)  yields 

[rJ(h),  f(h.)].  -b.(h.)/7jh).  (2) 

Hovever  {k.h.l,   2)  imply  that 

7*^  (h.)  -^   7    (h),  (3) 

hence 

-b.(h.)/7J(h.)  ~  [rJ(h),  f(h.)].  (^). 

By  using  (h)   and  (4.U.2,    3,   8)      in  (U.2.6)  we  obtain 

u(t,X)^(|A''/    zIexp|iX(V^   -ht)-Ti.(h)tj    [r^(h),f(h.)]rJ(h)dK.        (5) 


Here  we  sum  over  every  real  zero  ci)  =  h(K)  of  det  G  =  0,  then  over 

j   =  l,...,q_;   where  q  is  the  nullity  of  G(h),  (i.e.,  the  multiplicity  of  the 

eigenvalue  y{w)    of  G  which  vanishes  at  w  =  h). 

From  the  analog  of  (i|.  1,5), and  from  (4.1.13)  we  obtain 

t^ 
f(h.)  =  f(h.,K)  -  /  dY  rdT  exp  |-ix(k^y^-hT)+ n^(h)T  j  £  (t,Y)  .     (6) 

0 
Here  we  have  again  used  (i+.l|.2,  8  ).   We  have  also  used  the  fact  that,  by 

assumption,  f(t,X)  =  0  for  t  <  0.  We  now  insert  (6)  in  (5).   For  t  =  t^ 
this  yields 
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X  [rJ(h),  f(T,Y)]rJ(h).  (7) 

But  u(t  ,X)  must  be  independent  of  those  values  of  the  source  function 
f(t,X)  for  vhich  t  >  t  ^  hence  we  may  remove  the  "temporary  condition" 
(^.1.13).   Finally  since  t  was  an  arbitrary  positive  number,  we  obtain 

t 

X  \v^{K),   f(t,Y)l  rJ(K);    t  >  0.    (8) 

Here  we  sum  over  every  real  zero  co  =  h(K)  of  det  G  =  det[k  A  -  w(5'(w)]=0, 
then  over  j  =  l,...,q.   For  each  h(K) ,  q  is  the  nullity  of  G,  and  the 
eigenvectors  r'^(K)  =  r*^  [h(K)  ,k]  =  r  (h)  satisfy  the  equations 

G  rJ  =  0,   j=l,...,q;      G  =  k^A^'  -  h  g  (h)  (9) 


(r\  A°(h)  rJ)  =  &   •      i,j=l,...,q;  (lO) 


and 


(r^,  hB(h)  r^)  =  n  (K)6   ;      i,j=l,...,q.  (ll) 

The  last  equation  also  determines  the  values  of  t\  .(K)  =  ri  .  [h(K)  ,k]  =  tj  .(h) 
in  terms  of  the  matrix  Q  .        It  can  be  shown  that  if  the  numbers  tj  , .  . .  ,ri 
are  distinct,  equations (9 -11)  determine  the  eigenvectors  r  uniquely 
except  for  a  unitary  factor  expfia.]   (a.  real).   Such  a  unitary  factor 
has  no  effect  on  the  value  of  (8).   If  the  t)  .  are  not  distinct  the 
eigenvectors  are  uniquely  determined  by  (9-11)  except  for  a  imitary 
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transformation  (of  the  form  r  =  / a.  .  r'^ ,   where  (a.  .)  is  a  unitary  matrix) 

which  again  leaves  the  value  of  (8)  unchanged.   Since  (9),  (10),  (11) 
are  identical  to  (3-2.3  ,  ^,   H)  we  may,  for  the  case  of  constant  coefficients, 
identify  the  eigenvectors  of  chapter  3  and  chapter  k.      For  this  reason  we 
have  used  the  same  notation  in  both  chapters. 

Our  asymptotic  solution  (8)  of  the  problem  (U.1.1,  2,  ^4-)  is  the 
main  result  of  this  chapter.   It  was  derived  under  conditions 
(1-6)  which  are  listed  in  appendix  G.     For  real  source  functions  f , 
(8)  can  be  simplified  provided  we  assi^me  that  the  matrices  A  and  E  are 
real  (condition  9),  while  E   (w)  and  aj£)(co)  are  even  functions  (condition  10). 

We  recall  that 

G(a),K)  =  k^A^"  -  w£(co).  (12) 

From  condition  10  it  follows  that 

G(-c.,-K)  =  -G(c^,K),  (13) 

hence  det  G(w,K)  =  0  if  and  only  if  det  G(-(jJ,-_}^)  =  0.   From  this  we  see 
that  the  values  of  the  multiple -valued  function  co  =  h(K),  which  is  the 
solution  of  det  G  =  0,  can  be  so  chosen  that  for  each  value 

h(-K)  =  -h(K).  (14) 


We  now  define  the  fimctions  u,(t,X)  and  u  (t,X)  by  the  right  side  of  (8), 

except  that  for  u  we  sum  only  over  the  ■»    ^.   [  roots  co  =  h(K)  .   Then 
^±  I  negative  I  ^^ 

(assuming  that  we  may  neglect  contributions  from  any  zero  root)  we  have 
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u  f^  u.  +  u  .  (15) 

From  condition  9  '^^   see  that  we  may  choose  the  eigenvectors 
r*^  to  be  real,  and  since  G[h(-K),  -K]  =  -G[h(K),K],  they  may  be  so  chosen 
that 

r.(-K)  =  r.(K).  (16) 

From  (11),  (1^)  and  condition  10  we  see  that 

r,.(-K)  =Ti.(K).  (17) 

We  now  change  K  to  -K  in  the  integral  expression  for  u  ,  and  make  use  of 
(Ih) ,    (16),  and  (17).   The  result  is  ^ 

t    

u^it,X)    =i^      l^dYjdKJd^     )_       exp|ix[-k^(x^-y^)+h(K)(t-T)]-Ti^(h)(t-; 

j=l,. .. ,q 

X  [r'^,f(T,Y)]r'^;   t  >  0    (I8) 

But  from  {ih)   we  see  that  h(-K)  >  0  if  and  only  if  h(K)  <  0,   Therefore, 
by  taking  the  complex  conjugate  of  (I8)  we  find  that,  for  real  f 

u*(t,X)  =  u_(t,X).  (19) 

It  follows  from  (15)  that 

u(t,X)  ^  u,  +  u*  =  2  Re  u"*'(t,X)  .  (20) 

This  equation  holds  under  Qonditions9  and  10  provided  f  is  real,   u  (t,X) 
is  defined  by  the  right  side  of  (8),  except  that  the  sum  is  taken  only 
over  positive  roots  w  =  h(K) . 
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(5.1) 


5.   Radiation  from  soixrces 
In  this  chapter  we  consider  the  problem 


n 


v^  +X.a'',u^  =  f(t,X),    v(t,X)  =  /  F(T,X)u(t-T,X)dT;  (1) 

rt*.    v=l      V 


t    ^— r     '^X, .    ^   'a/v  '      ^   !v^     ,  /    ^      '^rn}  '/VL 


u(t,X)  =0     for  t  <  0.  (2) 

Here  f(t,X)  =   0  for  t  <  0  and 

The  matrix  t  and  the  hermitian  matrices  A  are  smooth  functions  of  X  in  all 
of  X-space,  and  t  satisfies  conditions  1-6.   The  source  function  f  may 
depend  also  on  X.   It  is  possible  to  obtain  the  leading  term  of  the 
asymptotic  expansion  for  X  ->  oo  of  the  solution  u(t,X)  of  (1,2)  for  a 
variety  of  types  of  source  functions  f .   In  order  to  illustrate  our  theory, 
and  in  particular  to  demonstrate  the  "indirect  method"  outlined  in  section 
3.11,  we  shall  consider  two  representative  examples  in  this  chapter.   Others 
are  discussed  in  \l,2,k,'^ ,16] . 

5.1  Rapidly  var.ylng  source 

The  source  function  f  will  be  called  "rapidly  varying"  if  it  is 
of  the  form 

f(t,X)  =  x^''^  g[x(t-T),  x(x-x^)J  (h) 

where  t  is  a  positive  constant  and  X  is  a  constant  vector.  An  important 
example  of  (h)    is  obtained  by  setting 

g(t,X)  =  5(t)6(x  )  ...  6(x  )g  =  &(t)6(X)g  (5) 

/*~      '~>-  -L  no  j,y^    rtJJ 


(5.1) 

where  g  is  any  constant  vector.   Then,  since  A5(AX)  =  5(x)  (k)   becomes 
X(t,X)  =6(t-T)6(X-X^)g^.  (6) 

In  general  we  shall  assume  that  the  function  g(t,X)  has  "compact 
support",  i.e.,  vanishes  outside  of  a  bounded  region  in  tjX-space.-^ 
Then  as  X  ^  oo  the  support  of  f(t,X)  shrinks  to  the  manifold  ''(consisting 


.^^^      /V^. 


of  the  single  point  (t,X)  =  (t,X  ).   ¥e  now  proceed  as  outlined  in  section 
3.3-   Since  ii(  consists  of  a  single  point,  s  has  some  constant  value  s 
on  l'^  .   This  constant  value  may  be  absorbed  into  the  undetermined  function 
z  therefore  it  is  convenient  to  take  s  =0.   Since  the  dimension  r  of 

w>  o 

T(|_is  zero,  there  are  no  conditions  (3'3'7)  and  we  may  choose  as  our  parameters 

r^  =  k^^;  J  =  l,...,n.  (7) 

This  yields  the  initial  conditions 


X  =  x^(r)  =  x^,   t  =  T(r)  =  -x,      K  =  K  (r)  =_r,  s  =  s^(r)  =  o.   (8) 


^  'w.     _  O  ^^  .,„    .P^/w     ,w'  O 


The  initial  value 

^  =  w  (r)  (9) 

O  /v\. 

is  determined  by  the  dispersion  relation  (3.2.1)  or  (3.2.2)  in  the  form 

det  G  =  det  \y  A^ {X   )  -  w  ?(lo  ,X  )1  =0  (lO) 

or 


(ll)  is,  in  general,  a  multiple -valued  function  of  the  vector _r_  which  may 


be  thought  of  as  a  point  in  n-dimenslonal  space ,  E  .   We  may  make  co  a 


single -valued  function  of  r,  by  letting  r  be  a  point  in  a  parameter  space' 


*This  assumption  can  be  weakened.   It  is  sufficient  to  assume  that  g 
decays  sufficiently  rapidly  at  infinity  for  the  integral  (3I)  to  exist. 


(5.1) 
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consisting  of  several  duplicates  of  E  ,  one  for  each  value  of  (ll).* 
By  this  simple  device  ve  conform  to  the  presentation  of  chapter  3,  in 
which  X  ,  T,  K  ,  s  ,  and  w  were  assumed  to  "be  single -valued  functions  of  T. 

/vS*      /wf"    o         o  z,^ 

From  (3- 3-1^)  s  is  given  by 

t 
s(t;r)  =  f   j?(t';r)dt'  (12) 


where 


and 


Jf(t;r)  =  k  (t;r0gjt;r)-u)  (r)  (13) 


^v^^'i:)  =\  l^^^'i:)' ^(*'i:)]5  v  =  i,...,n.  (lU) 

V 

Here  X(t;r),  K(t;r)  is  the  solution  of  the  ray  equations  (3.3-1)  with 

^^^     ^^w    ^^^     ^^^» 

initial  conditions  (8).   From  these  equations  we  find,  as  in  section  S-^j 
that 

^(^jr)  =  x^^  +  (t-T)hj^  (r,x^)+  o[(t-T)^,  (15) 


V    'V^ 


and 

.n+l-i 
^,x^;j  +  uut-T, 


j(t;£)  =  (t-T)"  det[\  ^  (£,X^)]  +  oHt-r)''''^].  (l?) 


It  follows  that  V  =  n  and 

T(T;r)  =detrh   (r,xjj.  (i8) 

We  recall  that,  in  chapter  3,  we  were  able  to  solve  the  transport 
equations  if  the  nullity  q  of  G  is  1  or  2.   If  q  =  1,  z_(t)  =  zj^t;!^)  is 
given  by  (3. 10. 7).    If  q  =  2  then  _z(t)   is 

*To  be  more  precise  each  space  E  is  restricted  to  those  values  of  P  for 
which  the  corresponding  value  of  (ll)  is  real.  "^^ 


(5.1) 
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given  by  (3. 10. 3,  k,   6)  provided  conditions  6,7,8,9  are  satisfied.   In 
either  case  it  is  necessary  to  determine  z(T;r)  vhich  appears  in  these 
formulas  for  z.   Let  us  suppose  that  we  have  determined  z.   Then  with  s 
given  by  (12),  the  asymptotic  solution  of  (l,  2,  k)    is  given  parametric ally 

by 

u~  e^^^^^'^C)  z(t;r);  (19) 


/v/>. 


X  =  X(t;r)  (20) 

(As  in  section  (3. 11),  for  each  point  (t,X),  (19)  is  to  be  summed  over 
all  values  of  r  which  satisfy  (20),  i.e.  over  all  rays  through  (t,}<^.) 

To  determine  z  we  use  the  "indirect  method"  outlined  in  section 
3.11.  We  first  consider  the  "canonical  problem"  with  constant  coefficients, 
A  =  A  (X  )  and  <t(w)  =  C^  (u),X  ).   For  this  problem,  we  find  (as  in  section 
3.6)  that  the  solution  of  the  ray  ecjuations  is  given  by 

\(^'X}  =  ^'  ^(^'£)  =  ^ov-^(*-^)gv'  ^v  =  \E5jp]'  ^21) 

and  (13)   becomes 

il(t;^)    =  7^g^    -  ^^,  (22) 

with  o)  given  by  (ll).   Now  (12)  takes  the  form 

■  s(t;r)    =   (t-T)SI  =  (t-T)[7^g^  -  coj.  (23) 

As  pointed  out  at  the  end  of  section  3. 10,  z  is  now  given  by  ($.10.9),  i.e. 


z(t) 


f[^]''exp(-(t-T),)  z;(t).  (21.) 


We  see  from  (21)  that  (I5-I8)  now  hold  exactly  with  no  remainder  terms, 
hence 
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.j(t)_ 


1/2 


=  (t-T) 


-n/2 


(25) 


Thus  (19)  and  (20)  become 


u~  exp  (i^[(t-T)(7^g^-w^)l-(t-T)Tl])  (t-T)""/^  z; 


X   =  X    +  (t-T)g  . 

V    ov   ^    '^v 


(26) 


(27) 


The  next  step  in  the  indirect  method  is  to  compare  (26,  27)  with 
the  asymptotic  expansion  of  the  exact  solution  of  the  same  problem.   To 
do  this,  we  first  insert  (l^■)  in  (4.5.8).  If  q.  =  2  we  apply  condition  7-  Then 

n-L  =  Ig  =  '^'  and 

0 


u-  r 


[^    I^f^^f^"   ^  exp{iX[7^(x^-y^)-h(t-a)]-n(t-a)} 


^^i^^^^}'    g^('^-T).  ^(l-'^o^l}^(^'  *  >  0-  (28) 

Here  we  sum  over  every  real  zero  w  =  h(r;X  )  of  det^y  A  (X  )-wc.('jJ,X  )^  =  0, 
and  then  over  j=l,2ifq  =  2.   Ifq  =  lwe  may  omit  the  index  j.   ¥e  now 


introduce  the  transformation 


,  n. 


t'    =  Ha-y)  ,     X'    =  X(Y-X^),   dt'    =  Xda,      dX'    =  r  dY, 
in  (28).      The  result  is 


/v2 


(29) 


u 


-{^     /dr2Lexp   {ix[7,(x^-x,,)-h(t-T)]-(t-T),)X[rJ(r),a,(h,r)]rJ(r) 

*      '     '■'  u=h  j 


where 


(30) 
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(5.1) 


X(t-T) 

^(h,r)  ^Jdx^l         dt'  exp  {i[ht' -  7^x;]-hit'/x}g(t' ,x:) 

-Xt 


0         for  t  <  T 


< 


I  dX-  J     df    exp  {i[hf  -  7^x;])g(f  ,X:)      for  T  <  t 

-00 

We  note  that  (the  asymptotic  form  of  )  a  is  simply  the  foiorier  transform 
of  g.  It  is  independent  of  X,  For  the  special  case  (6)  we  see  from  (5) 
and  (31)  that  a(h,r)  —   g  . 

We  now  apply  the  method  of  stationary  phase*  to  obtain  the 
asymptotic  expansion  of  the  integral  (30):  We  introduce  the  function 


(31) 


h(r)  =  y   (x  -X  )-h(t-T) 


Then 


i>       =  X  -X  -(t-T)g  ; 

'  V 


V   -^ 


(32) 


(33) 


Stationary  points  are  determined  "by  the  condition  that  (33)  vanish  for 
V  =  l,...,n.   At  the  stationary  points, 

\  =   ^ov  +  (^-'^)gv^  (3^+) 


^    =    (t-T)[7^g^-^J, 


w  =  h(r;X  ) , 


and 


V  7 


V    n 


(35) 


(36) 


Therefore,  for  t  <  t  the  stationary  phase  formula  yields 


*A  simple  discussion  of  this  method,  in  n  dimensions,  is  given  in  (^10 ^ 
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u-.(i^"^.-.,-"/^s:[i-(.,,;)i]-/^ 


X  exp{±x[{t^)fy^g^   _uj^)]_(t-T)  -±n/k   sig(h^  k  )}  X  k^fL(h.i:)]r"^-   (3?) 

Here  sig  (h,  .  )  denotes  the  signature  of  the  matrix  (h,  ,  ).   At  each 
^  k  k  ^  k  k 

point  {t,X)   we  must  sum  over  all  r  satisfying  (3^)'   Alternatively  if  we 
adopt  the  extended  definition  of  the  parameter  space  introduced  earlier  in  this 
section,  we  may  omit  the  first  summation  sign  in  (37)-   It  is  then 
understood  that  (37)  is  to  be  summed  over  all  T   in  (j  that  satisfy  (3^)- 
We  now  compare  (3^,  37)  with  (26,  2?).   ¥e  find  that  they  are 
identical  provided 

n/2  "^"/^  sig(hj^  ^   )    q_ 

3x;r)=(|^)   [|det(h,,)|  -V^e  ^^   Z  k^a(h,r)lrj.   (38) 

Here  q  =  1  or  q  =  2,  h  =  h(r,X^)  and  h^^  ^     =  ^^-^  h(K,X^)|      -where 

oj  =  h(K,X  )  is  the  (multiple -valued)  solution  of  det[k  A  (X  )-toC.(co,X  )J  =  0. 
Furthermore,  we  see  from  (ij-. 5.9-11)  that  the  vectors  r  =  r  (P)  satify 


v\-        ^^/"^        rv^ 


the  equations 

Gr-^  =  0;    j  =  l,...,q;    G  =  y  ^  (X^) -h  ^{h  ,X^)  ;  (39) 


and 


(r^,  A°(h)r-^)  =  S   ;      i,j=l,...,q;  (kO) 

[r_^,  hc?)(h,X^)r'^]  =  Ti(r)5.  ;        i,j=l,...,q.        (^l) 


These  equations  determine  the  r  except  for  a  unitary  transformation  which 
leaves  (38)  invariant. 
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¥ith  z'  given  by  (38)  we  may  ncfw,  in  principle,  complete  the 
solution  of  (1,  2.,    h) .      Of  course  this  requires  the  solution  of  the  ray 
equations  which  is  not,  in  general,  possible  by  analytic  means.   Under 
special  assumptions  about  the  functions  A  (X)  and  C^  (f^,X)  this  is  sometimes 
possible*  and  our  method  leads  to  escplicit  formulas  for  the  asymptotic 
solutions.  ¥e  shall  not  carry  out  the  details  of  this  procedure  here 
for  any  problem  with  non-constant  coefficients. 


*See,  e.g.,  [l  J 
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5.2  Oscillatory  source 

In  this  section  we  shall,  for  convenience,  restrict  our  attention 
to  three-dimensional  space  (n  =  3)  and  we  shall  assume  that  the  dispersion 
relation  is  isotropic ,  i.e. 

(J  =  h(K,X)  =  h(k,X);        k  =  Ixl .  (l) 

We  further  assume  that  the  multiple -valued  function  (l)  has  (for  each  X) 
a  single -valued  inverse 

k  =  m(a))  =  m(oj;X)  ;       m' (oj)  /  0.  (2) 

The  function  (2)  is  defined  for  w  in  some  subset  Q(X)    of  the  real  axis. 

•ws. 

For  each  X,  f2(X)  consists  of  all  (real)  values  which  (l)  assumes  when  K 
takes  on  all  real  values  (or  k  takes  on  all  non-negative  real  values). 
We  shall  see  in  chapter  7  that  the  above  assumptions  are  suitable  for 
electromagnetic  waves  in  isotropic  media. 

We  now  consider  source  functions  of  the  form 

f(t,X)  =  x3g[x(X-Xp),t]  e-^^%^^\  g  a  0  for  t  <  0.     (3) 

If  we  take 

g(X,t)  =  5(X)  g  ft)  (h) 

then  (3)  becomes 

f(t,X)  =  6(X-X  )g  ft)  e-^^^o(t)  ■  (5) 

*«   /*\      /«  ^o   o  • 

(5)  represents  an  oscillatory  point  source.   If  g  is  constant  and  q^(t)=  ilA, , 

the  source  is  time-hamionic  with  frequency  ^X.  If  g  is  not  constant  but 

a  (t) =  •^  we  may  call  it  an  amplitude -modulated  time-harmonic  source  of 

frequency  ^X. 
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We  assime  that,  in  general,  g(X,t)  has  compact  support*  for  each 
t.   Then  as  X  ^oo  the  support  of  (3)  shrinks  to  the  one -dimensional  manifold 
^Y  ,   defined  parametrically  (with  parameter  t)  ty 

X  =  X  ,   t  =  T  ;         T  >  0.  (6) 

-VIA     /W.O  — 

From   (3)    it   is  reasonable   to  assxjme  that  only,    s   is  given  by 

^[>.XJ=  -qjT).  (7) 

We  shall  verify  this  assumption  shortly  by  comparison  with  the  solution  of 
the  canonical  problem.   Differention  of  (7)  yields 

0)  =  4(t).  (8) 

Proceeding  as  in  section  3-3  we  choose,  as  independent  parameters,  T  =   {x,ci,y), 
where  a  and  7  are  polar  angles  (latitude  and  longitude)  defined  by 

fcos  a,  sin  a  cos  7,  sin  a  sin  7J  =  A  =  K  /k  ;   k  =  |K  1.   (9) 

The  initial  values  of  X,t,s  cj,K  are  then  given  as  functions  of  F  by  (6), 
(7),  (8),  and 

K  =kA;        k  =  m(co;X  )  =mrq  (t)  ;X  "I  .  (10) 

Since  m(cj)  is  single-valued,  all  the  initial  values  are  given  as  single- 
valued  functions  of  F,  as  required  in  section  3'3j  and  the  parameter  spaceyj 
consists  of  all  values  of  F  for  which  4Jt)  is  in  n(X  ),  t  >  0, 
and  the  polar  angles  are  in  the  intervals  0<a<n,  C)<7<  Sjt. 

We  now  denote  by  X_  =  X(t;r),  K  =  K(t;r)  the  solution  of  the  ray 
equations  (3. 3.1).   (For  each  value  of  F  =  (T^a,7)  we  use  in  the  ray 
equations,  that  value  of  the  multiple -valued  function  (1)  for  which 
hrm(c3  ,X  );Xl=q  (x).)  We  note  that  since  oj  is  constant  (8)  holdsalong  every 
ray,  and  we  set 


*It  is  sufficent  to  assume  that  g  ^  0,  as  |x|  ->  00,  sufficiently  rapidly  so 
that  the  integral  (30)  converges. 
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A(t;r)  =  (a^,  3.^,   a^)  =  k"Vt;£)  =  (m[oj;X(t;r)])"^  U^'>'Q  >  (H) 

g^  =  h   =  h'(k)a^  =  [m'(a))]-\^;  G  =  (g  ,g  ,g  )  =  (m'  ['^;X(t;r)]}  "\(t;r) , 
V  J  -^ 

and  (12) 

m[w;X(t;r)l 
je(t;r)  =  K-G  -  0.  =        -^         -  ^.  (13) 

m-  [c^;x(t;r)] 


0^  /^'v-  /fc^ 


Then 

t 

s(tjj^)  =  -q^(T)  +  /'j|.(t';r)dt'.  {Ik) 

T 

As  in  section  3-^  we  see  that  the  ray  equations  may  be  used  to 
obtain  the  expansion 


X(t;r)  =  (x^,x^,x_)  =  X  +  — [cos  a, sin  a  cos  y,    sin  a  sin  y\ 

-    -        1   2   3      _o    ^.[4j^).y 

+  o[(t-T)^].       (15) 

By  differentiating  (15)  we  find  that* 

ax   f^^      Sxg   atxN      [l+(t-T)m"q^/m'J 

^  =  1-5^,  -5p,.^)  = -; L^os  a,  sin  a  cos  7,  sin  a  sin  y] 

+0(t-T),  (16) 

SX       /Sx        ctx^      ^  \  2 

^  =  I  a^'  a^'  aTJ  =  ;^t-"^^  ^^  "°"  ^  "°" '''  "°'  ^  '^^  7]-K)[(t-T)  ], 

(17) 

and 

SX       /cDC-|^     Stx^        ctx  \      ^  2 

^  =  1^,  -^  ,   •5r-)=  ^  [Oj  -sin  a  sin  y,    sin  a  cos  7]+  o[(t-T)  ].    (I8) 

From  these  equations  it  is  easy  to  compute  the  jacobian. 


*See  footnote  on  the  next  page. 
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/^l 

^o 

St 

1 

^  ^ 

^ 

j(t;r)  =  det 

^1 

^2 

^ 

\ 

V 

^2 

(19) 


The  result  is* 


j(t;r) 

/V\. 

— 

(m 

^  (t. 
•)3 

or 

J(t, 

;r) 

=  - 

-  sin  a 

It 

follows  that  1 

^  = 

2  and 

j(t 

;r) 

=  • 

-sin  a 
(in')3  ' 

(t-T)^[l+(t-T)m"(^/m']  +0[(t-T)3l  (20) 


{t-xf   +  0[(t-T)3].  (21) 


-  sm  a 


(m'[4^(T);xJ)- 


(22) 


Just  as  in  section  5-1  we  can  now,  in  principle,  complete  the  asymptotic 

solution  of  our  problem  if  q  =  l,or  if  q  =  2  and  conditions  6,7,8,9  are  satisfied; 

except  that  z  =  z(T;r)  is  undetermined. 

"■^       ^w%.       -v\. 

In  order  to  determine  z,  we  again  consider  the  canonical  problem. 
We  set  A  =  A  (X  )  and  ^  (to)  =C(w,X  ).   Then  m(io)  =  m(u);X  ),  and  on 

-v^V.  -\.«^  __,.^ 

each  ray. 


CJ: 


4   (t),   k  =  m(oj)  =  m  [6    (t)  I  ,   K_  =  kA,  A_  =  A     =    [cos  a,    sin  a  cos  7,    sin  a  sin  7] 

(23) 


*In  (16)  and  (20)  the  term  (t-T)m"q/m'  could  be  omitted  since  it  is  of  the 
same  order  as  the  remainder.   We  keep  it  to  facilitate  our  later  discussion 
of  the  problem  with  constant  coefficients. 
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Furthermore 

G^=A/m'(a3),    s  =  -q^(T)+(t-T)£  ;      ^  =  [m(u))/m' (w)] -u,       {2h) 

and  the  rays  are  given  by 

X  =  X  +  (t-T)A/m'(oj).  (25) 

We  now  find  that  (l6-20)  are  valid  for  all  t  >  t  with  the  remainder 
terms  equal  to  zero.   Hence 

1(t) 


_j(t) 
Wow  (3.10.9)  yields 


=  (t-T)-^[l  +  {t-x)m"Cq_JiJm'ir^^)\~^^^.  (26) 


^(t)  =  (t-T)"^[l+(t-T)m"(q^ym'(q^)|^'^~exp|-(t-T)Vl  z(t)      (2?) 
and^  from  (3-11 .1^2)  we  have 

u   ^  (t-T)-^[l+(t-T)nf(q^)^Mq^)  r^/2  exp  ^^[(t-T)  (^U>qo(T)] -(t-T)Tijz(T) , 

X  =  X  +  (t-T)A/m'(qQ).  (28) 

We  must  now  compare  this  result  with  the  asymptotic  expansion 
of  the  exact  solution  of  the  canonical  problem.   Therefore  we  insert  (3) 
in  {h.3.8).     We  then  introduce  the  transformation  X'  =  X(Y-X  ) ,  dX'  =  X^dY, 
of  the  integration  variables.   Thus  we  obtain 

0  ^ 

X  [r'^(K),  a(K,T)Jr'3(K),  (29) 

where 

-ik  X' 


/-±ts.   X 
e   "^  ""  g(X',T)  dX'.  (30) 
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(5.2) 


In  (29)  we  sura  over  all  values  of  (l).   If  q  =  Iwe  omit  the  index  J .  If  q 
we  Slim  over  j  =  1,2.   In  the  special  case  {k)   we  find  that  a(K,T)  =  g  (t). 
We  now  apply  the  method  of  stationary  phase  ["  lo]  to  the  U-f old  integral 
(29):  We  introduce  the  function 


=  2 


=  k  (x  -X  )  -  h(t-T)  -  q  (xj 
V^  V   vo  o 


stationary  points  are  determined  by  the  conditions 


(31) 


k    =  i\-\o^-i^-r)&y  =  0;^=^,2,3;      U^  =  \    =h'(k)a^ 


k^k);   (32) 


and 


t)^  =  h-q^(T)  =  0. 


At  the  stationary  points, 

(j)  =  (t-T)[k^g^-h]-q^  (t)  =  (t-T)[kJi'(k)-h]-qj;T) 
and  the  second  derivatives  of  i  are  given  by 


(33) 


(3^) 


k        =  h,      =  h'  (k)a   ; 
^^k  T  k  ^    '    v' 

V  V 


t„^-3„(T); 


(35) 


and 


k  k     =   -^^-^K  k     =   -(t-T)[h"(k)a^a     -^  ^  i\^-,-^)]'  06) 


We  now  introduce  the  matrix 

Vi  ^^1^2  V3 

^^2^1  V2  ^^2^ 

K^k^        't'k^k^  i3k3        i^T 


^Tk, 


^. 


rK 


^k. 


(37) 


'TT 
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By   choosing   a  co-ordinate  system  in  which  A  =  (a,  ,  a  ,  a  )  =  (1,0,0) 

r>~     J.    ^    3 

(35)  and  (36)  simplify  greatly,  and  it  becomes  easy  to  evaluate  the 
determinant  of  ( 37 ) •   We  find  that 


.2,,  .,h 

(h-) 


aetO  =  -fety^  £l-(t-T)h^qn].  (33) 


In  order  to  determine  the  signatiore  of  (o)  we  may  find  its  eigenvalues. 
They  are  the  solutions,  x,  of  the  equation  det($  -  xl)  =0.   This  deter- 
minant is  also  easily  evaluated,  and  the  equation  for  x  factors  into  the 
two  equations 


(t-x)  I'  +  X 


2 

=  0,  (39) 


and 


[(t-T)h"    +  xj  [q^+x]  -    (h-  f    =   0.  (i^O) 

The  two  roots  x  =   -(t-T)h'/k  of   (39)   vanish  at  t   =  t,   and  for  t  >  t 
sgn  x  =   -sgn  h'    =   -sgn  m'[q  (t)].      The  two  roots  of   (i|0)    do  not  vanish 
at  t   =  T  and  are  given  by 


2x  =   -^±    y(q/+  Mh')2    .  {hi) 

It  is  clear  that  one  of  these  roots  is  positive  and  the  other  is  negative. 
Hence  in  a  neighborhood  of  t  =  t  (t  >  t);   i.e.,  up  to  the  next  caustic 
point,  sig  $  =  -2  sgn  h'  =  -2  sgn  m'(q.  ). 

The  stationary  phase  formula  now  yields  for  (29) 


^^^z{(t-.^^V-^^]} 


-1/2 

X  exp|i\[(t-T)(kh'-h)-qJ-Ti(t-T)-in/2  sgn  h' I  IV<^(K)  ,a(K,T)]rJ(K)  , 
h=  4o(t),  x^  =  x^^  +  (t-T)h'(k)a^;  J£  =  ^  =  1^(3^,32,3^).      (42) 
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Since  h'(k)  =   , /  v   and  h"(k)  =  -  — ^— ^  ,   this  becomes 
m'(c.)  (^,^3 


>   X-       1  o  ^""  -1/2 


1  +  (t-T)  ^o 

m' 


exp|ix[(t-T)(^,-q^)-qJ-(t-T)Ti-irt/2   sgn  m'j  [t^^Ck)  ,a.(K,T)J  r^(K) ; 

K=m(4o)A;   X  =  X  +  (t-T)A/m' (q  ) .  (43) 

In  (U2)  u(t,X)  is  represented  parametrically  with  parameters  t,   and 

K  =  (k-,  ,  k^,  k  ).   These  parameters  are  not  independent,  because  they 
■•^  123 

are  related  by  the  equation  h  =  qQ(T).   In  {k2)    the  indicated  sum  is  to  be 
taken  over  the  values  of  (l)  as  well  as  over  the  index  j   (if  q  =  2). 
In  (^3),  u(t,X)  is  represented  parametrically  by  3  independent  parameters. 
These  pareameters  are  t  and  two  angles  which  determine  the  unit  vector  A. 
They  may  be  taken  as  in  (23).   In  (^3)  the  indicated  sum  is  to  be  taken 
only  over  the  index  j  (if  q  =  2). 

We  may  now  compare  (28)  and  (U3).   We  see  that  they  agree  exactly 
if  we  set 

z(T;r)  =  ^  m(4,)k(q,)]2e-'^/2  ^^   ^'^'^^  {r^  (K)  MK,.)]rJ  (K)  .  (kk) 

J=l 

Here  K -m[qQ(T)jA^,   A  =  fcos  a,  sin  a  cos  j,    sin  a  sin  7},   and  T  =   {x,cl,7)' 
As  in  section  5-1,  (44)  may  now  be  used  to  complete  the  asymptotic  solution 
of  the  problem  (5.1.1,2),  (5.2.3)  for  the  case  of  non-constant  coefficients. 
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It  is  interesting  to  note  that^  for  the  case  of  constant 
coefficients  the  jacobian  has  not  only  the  zero  of  order  2  at  t  =  t 
but  also  a  zero  of  order  1  at  the  caustic  point  t  =  t  >  r,  determined 
by 

Such  a  point  exists  if  and  only  if  h"q  >0.   Then  for  t  >  t  the  signature 
of  $  changes.   To  find  sig  $  we  write  (UO)  in  the  form 

x^  +  Bx  +  C  =  0;   B  =  (t-T)h"  +  q^,   C  =  (t-T)h"^-(h' )^.  (^6) 

Then 

2x 


-B  ±  /b^  -  ^C  ,  (hi) 


and 


sgn  C  =  sgn  (t-t  ).   Thus  the  positive  solution  of  (k'j)   becomes 


negative  and,  for  t  >  t  ,  sig  $  decreases  by  2.   It  follows  from  the 
stationary  phase-formula  that  {h3)    then  must  be  multiplied  by  a  factor 
e~  '  .   This  confirms  the  "phase-shift  rule"  of  section  3' 10  for  this 
case. 
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6.   Dispersive  and  Won-disperslve  Systems 

6.1  Non -dispersive  systems 

We  shall  say  that  our   general  system  of  equations  is  non- 
dlspersive  if   <5=  A(X)  is  independent  of  co.   If,  furthermore, 
O  =  C  -(ix)  J)    where, 

£(co,X)  =A(X),     and       '^(^.X)  =  C(X)/w,  (l) 

then,  as  in  section  2.2,  we  see  that  the  general  system  of  equations 
(2.1.26,  27)  becomes  the  symmetric-hyperbolic  system 
n 

A  Ut  ■*■  ^  ^  ^x  "^  ^  ii  =  ^-  (2) 

v=l      ^V 

We  also  see  that  A  =  (^8),^  =   A  and  the  dispersion  relation  becomes 

det  G  =  det  (k  A^  -  wA)  =  0.  (3) 

By  condition  2,  A  must  be  positive -definite. 

It  is  now  easy  to  verify  that  the  dispersion  equation  (i.e. 

the  partial  differential  equation  (3)  for  s,  with  co  =  -s^  and  k  =  s   ) 

t       V     X  ' 
V 

is  the  characteristic  equation*,  the  level  surfaces  of  s(t,X)  are 

characteristic  hyper surface s»,  and  the  rays  are  the  bicharacteristics*  of 

(2). 

If  we  express  the  dispersion  relation  (3)  in  the  equivalent  form 
^  =  h(K;X)  (4) 

we  see  from  (3)  that  the  (multiple -valued)  function  h  is  homogeneous  of 

degree  one  in  K:  i.e.  for  all  real  e 

*See  [  3  ]. 
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h(eK;X)  =  eli(K;X).  (5) 

From  the  Euler  equation  for  homogeneous  functions  it  follows  that 

/  =  k  h,   -  h  =  K-G  -  to  =  0,  (6) 

V 

hence  from  (3-3>3)  ^we  see  that  s  is  constant  on  bicharacteristics. 
This  verifies  the  well-knovm  fact  that  bicharacteristics  lie  on 
characteristic  hyper surfaces.   In  section  3-3  we  pointed  out  that  the  term 
"ray"  is  sometimes  used  for  the  projections  X  =:  X(tjr)  of  the  space-time 
curves  (t,X)  =  [jt>X('t;r)j  into  X-space.   For  non-dispersive  systems  we 
shall  refer  to  the  space-time  curves  as  bicharacteristics  and  to  their 
projections  as  rays.   By  differentiating  (6)  with  respect  to  k  we  find 

H" 

that 

n 


/   k  h,  ,   =0,  (7) 

v=l      V  ^ 

and  from  this  equation  it  follows  that  the  matrix  (h  ,  )  is  singular, 

k  k 

i.e. 

det(h  J.  )  =  0.  (8) 

V  \i 

Since  a  non -dispersive  system  is  a  special  case  of  our  general 

system,  the  results  of  this  paper  apply  with  very  few  exceptions,  and  we 

may  construct  asymptotic  solutions  of  (2)  according  to  the  general  theory 

of  the  preceeding  chapters.   Of  course,  the  parameter  1  no  longer  appears 

in  the  system  of  equations  (2),  but  does  appear  in  the  data  of  the  problem 

(i.e.  the  function  f,  initial  data,  boundary  data,  etc.).   The  asymptotic 

theory  of  (2)  with  f  =  0,  and  with  oscillatory  initial  data  of  the  form 

iXs  (X) 
^(o.x)  =^   °"  l°(x),  (9) 
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was  first  discussed  by  P,D.  Lax  [_  9  J  • 

It  is  interesting  to  note  that  although  the  non-dispersive 
system     is  certainly  simpler  than  the  general  system,  singularities  • 
appear  in  the  special  case  which  do  not  occur  in  general.   Thus,  for 
example  the  discussion  of  section  5-1  is  not  valid  for  non-dispersive 
systems,  for  then  we  see  from  (8)  and  (5.1.38)  that  ^  is  infinite.   On 
the  other  hand,  the  discussion  of  section  5-2  is  valid  for  non-dispersive 
(isotropic)  systems,  and  the  results  in  this  case  simplify  considerably 
by  virtue  of  the  fact  that  h"(k)  =  0. 

For  non-dispersive  systems  the  factor  j(t  )/j(t)  which  appears 
in  all  of  our  formulas  for  z  has  a  special  geometrical  interpretation 
which  is  most  easily  examined  in  the  case  of  3  space  dimensions  (n  =  3) : 
Since  s  is  constant  on  bicharacteristics, 

s[t,x(t;r)]  =  sjr).  (lo) 

Let  us  choose  otir  parameters   P  =   (7-,  ,7r^,7-~,)    so  that 

"^      J-   c:   j 

s,(r)  =3^(7^).  (11) 

Then  for  each  fixed  value  of  y    ,      (t,X)  =  rt,X(t;7  ,7  ,7  )1  is  the 
parameteric  equation  of  a  characteristic  hypersurface  with  parameters 
(t,7  ,7  )  and  for  each  fixed  value  of  (t,7  ),  X  =  X(t;7  ,7  ,7  )  is  the 
parameteric  equation  of  a  wave-front*  with  parameters  {yr^,y   ).     The 
vectors 

^1  =  #^      ^2-  ?2  ^''^ 

are  tangent  to  the  wave-front.   Now  differentiation  of  (10,  11)  yields 

K-X,  =  s'(7J,        K-X„  =  K-X-  =  0,  (13) 

/f  '"l    o  1  '        ,^  ^2       '"  "^i 


■^The  terra  "wave-front"   is  defined  in   section  6.3. 
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hence  K  is  orthogonal  to  X^  and  X^i  i.e. 

^2  "^  ^3  =  ^S^  (1^) 

where  b  is  a  scalar. 

For  fixed  7,  ,   let  7  and  7  vary  in  infinitesimal  intervals 
of  lenth  d7  and  d7  .   Then^  for  fixed  t,  the  point  X(t;7-,  ,7  ,7  ) 
varies  in  an  infinitesimal  region  on  the  wave -front  of  area  da,  and 
the  rays  vary  in  an  infinitesimal  "tube"  of  rays  of  cross-sectional 
area  da.   Since  the  vector  G  is  parallel  to  the  rays  it  is  easy  to  see  that 

and 

^^  =  2^2  ^  -53-K  -f-^  =  bkd72d73.  (16) 

Here  g  =  |g|  and  k  =  |k|.   Wow,  by  definition, 

j(t)  =  j(t;r)  =  x^.x^  X  x^  =  b  K-x^  =  b  s;(7^).        (17) 

Hence  (6),  (15),  (I6)  and  (I7)  yield 

■(f\   _  lof   da    _  fo   da  ,  nv 

'^^    '  w  d.y^d.7       ~   k  d72'i73  '  ^      ' 

Since  La,s'  ,  ^.7  ,  and  d7  are  independent  of  t  we  find  that 


j(t^)   g(t^)   da(t^)   k(t)  da(t^) 


j(t)     g(t)   da(t)  -  k(t^)  da(t) 


(19) 


o 

Here  j(t)  =  j(t;r),  k(t)  =  k(t;r)  and  g(t)  =  g(t;r).   da(t)  is  the  area 
of  the  portion  of  the  wave-front  Intersected  by  the  tube  of  cross-sectional 


(6.2) 


area  da(t).   It  is  important  to  note  that,  since  7,  was  held  fixed, 

the  tube  of  rays  is  formed  from  rays  associated  with  a  single  character- 


istic hypersurface ,  i.e.,  projections  of  bicharacteristics  which  lie  on 
a  single  characteristic  hypersurface.  Since  G  is  parallel  to  the  rays  and 
K  is  normal  to  the  wave -front  we  see  that  if  0  is  the  angle  between  these 

/Wl 

two  vectors,  then 

da       ^   KiG    cj  ,      . 

r—   =    cos  ©  =  f*-^  =  :; (20) 

da  k  g   k  g  ^      ' 

in  agreement  with  (I8). 

For  the  non -dispersive  system,  either  of  the  expressions  (19) 

may  be  substituted  for  j(t  )/j(t)  wherever  it  appears  in  our  earlier 

formulas. 

6.2  Progressing  wave  solutions  of  the  non-dispersive  system 

We  have  seen,  in  section  3-3  that  the  (multiple -valued)  solution 
s(t,X)  of  the  dispersion  equation  is  determined  by  its  initial  values 
s  (r)  on  the  manifold  ''C.        Similarly  z  is  determined  by  z(T;r)  [or 
z^(T;r)  if  f^   is  a  caustic  |  .   Once  s  and  z^  are  found,  the  function 
exp{iXs}z  is  an  asymptotic  solution  of  the  general  system  of  equations. 

In  the  non-dispersive  case  it  is  easy  to  see  that  if 
exp{iXs}z,  X  =  X(t;r)  is  an  asymptotic  solution,  and  e  is  any  real 
number  then  the  function  (defined  parameter! ca3J.y  by) 

exp{iXes]3,      X  =  X(t;r),  (l) 

is  also  an  asymptotic  solution.   To  see  this  we  simply  multiply  s  ( r) 
by  e.   Then  since  s  is  constant  on  bicharacteristics,  s,  k  =  s   ,  and 
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to  =  -s,  are  changed  by  the  same  factor,  while  g  =  9(jj/Sk  and  the  bi- 
characteristics  are  unchanged.   Since  j(t  )/j(t)  is  unchanged,  z  is 
unchanged.  Now  if  e(x)  is  an  arbitrary  function,  e(e)  is  a  constant, 
hence 

?(e)  exp{i\es}z,       X  =  X(t;r),  (2) 

is  also  an  asymptotic  solution.   By  integrating  (2)  with  respect  to  e 
we  obtain  the  asymptotic  solution 

u  ~  erxs(t;r)]z(t;r),   "  X  =  X(t;r),  (3) 

which  involves  the  arbitrary  "wave-form"  function 


e(T)  =   /   exp{ieT}  e(e)  de.  (4) 

—  00 

(According  to  our  usual  convention,  at  each  point  (t,X),  (3)  is  to  be 
summed  over  all  values  of  r  for  which  X(t:r)  =  X,  i.e.  ,  over  all  bi- 
characteristics  which  pass  through  the  point.) 

Asymptotic  solutions  o'  the  form  (3)  of  non-dispersive  systems 
are  called  progressing  waves.   Applications  of  progressing  waves  are 
discussed  in  [ll,13J.  Here  we  mention  only  that  in  using  progressing 
waves  to  obtain  asymptotic  solutions  of  problems  for  non-dispersive 
systems  of  equations,  the  data  of  the  problem  yield  not  only  the  values  of 
s  and  z  (or  g)  on  /TZ  ,  but  also  the  wave -form  function  e(T). 
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6.3  Group  velocity,  phase  speed,  and  disturbance  speed  for  dispersive 
and  non-dispersive  systems 

In  this  section  we  return  to  the  consideration  of  dispersive 
systems,  but  we  shall  make  some  comparisons  with  the  non-dispersive  case. 
Many  dispersive  systems  satisfy  condition  11  of  appendix  G.   The  system 
of  electromagnetic  field  equations,  with  a  dielectric  permeability 
function  of  the  kind  discussed  in  appendix  D  is  of  this  type.   For  a 
system  which  satisfies  condition  11,  it  is  convenient  to  refer  to  the 
system  with  G  replaced  by  its  limiting  value  A  as  the  corresponding 
non-dispersive  system.   (The  latter  system  is  given  by  (6.1.2)  with  C  =  0. ) 
From  the  ray  equations  (3.3.I)  we  see  that  space  point  X  =  X(t;r) 
which  remains  on  a  given  ray  as  t  increases  moves  with  a  velocity 
X  =  G  =  (g-,,...,g  ).   The  vector  G  is  often  called  the  group  velocity 


(or  ray  velocity)  vector,  and  its  magnitude  g  is  called  the  group  speed 

(or  ray  speed) . 

From  (3.1.6)  we  see  that  the  vector  K  =  (k-,,...,k  )  is  identical 

to  the  vector  Vs  =  (s   ,...,s   ).   Hence  for  each  fixed  t  the  unit  vector 

1      n 

A  =  (a^,...,aj,     a^  =  Vk,      k^  .  K-K;  (l) 

is  normal  to  the  level  surfaces  of  the  phase  function  s(t,X).   These  level 
surfaces,  which  are  defined  by  the  equation 

s(t,X)  =  constant,   t  =  constant  (2) 

are  called  wave -fronts  (or  phase -fronts) .   As  t  increases  the  wave-fronts 
move  through  space.   Let  X  =  X^(t)  =  (x  ,  ...,x  )  be  a  point  which  is 
constrained  to  remain  on  a  wave-front  (2).   By  differentiating  (2)  with 
respect  to  t  we  find  that 


(6.3) 
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-w  +  ka  X  =  -w  +  k  X  =0.  (3) 

The  quantity 

p  =  X-A  =  X  a,  =  :;—  ,  (k) 

which  is  the  component  of  the  velocity  X^  of  the  constrained  point,  in 
the  direction  of  the  normal  to  the  wave -front,  is  called  the  phase -speed. 


From  the  dispersion  relation  (3. 2.1)  we  see  that  p  is  a  solution  of  the 
equation 

detCa^A^"  -  p<f )  =0.  (5) 

There  is  one  phase  speed  corresponding  to  each  value  of  h,  i.e.,  correspond- 
ing to  each  value  of  the  multiple -valued  function  s.   For  the  correspond- 
ing non-dispersive  system  the  phase  speeds  are  the  normal  speeds*  v  which 
satisfy 

det(a  A^  -  V  A)  =0,  (6) 

and  from  (6.1.6)  we  see  that 

A-G  =  p  =  V.  (7) 

This  equation,  which  asserts  that  the  component  of  the  group  velocity 
in  the  direction  of  the  nonnal  to  the  wave -front  is  equal  to  the  phase - 
speed,  holds  only  in  the  non-dispersive  case.   For  the  general  case  (subject 
to  condition  11)  it  is  replaced  by  an  inequality  which  we  now  derive: 

Since  a  A  is  hermitian  and  A  is  positive  definite,  there  exist 
vectors  s  ,...,s"^  such  that 

(a  a''  -  V  A)s"^  =  0;     (s^,  As"^)  =  6  ;      i,j  =  l,...,m.   (8) 


*  See  [3  J 
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V  ,■•  ..,v  are  the  roots  of  (6).   An  arbitrary  vector  r  may  be  expressed 
In  *^ 

in  the  form 


m 


P  .  s" 


It  follows  easily  from  (8)  and  (9)  that 


m 


m 


(r,  a  A^  r)  =  / v. p.,    (r,Ar)  =  /   p.. 


Hence 


(r,a  A^r)  <  v    (r.Ar);    v    =  max  v.. 

J 


^f^ '     y      ^ 


max: 


(9) 


(10) 


(11) 


Wow  let  r  be  one  of  the  null  eigenvectors  of  the  matrix  G  (section  3-2) 
From  the  basic  identity  (3.2.8)  we  see  that 


A.G  =  aR^  =   {r,   o. J^  v)  ,  (12) 

while  condition  11  and  the  basic  identity  (with  v  =  O)  imply  that 

(r,  Ar)  <  (r,  A°r)  =  1.  (13) 

(11),  (12)  and  (13)  now  yield  the  inequality 

A-G  <  V   ,  (1^) 

,^  '^  —    max 

which  is  the  required  generalization  of  (7)-   It  is  well  known  that  v    is 

max 

the  normal  speed  of  a  disturbance  front,  i.e.,  a  moving  surface  which  re- 
presents the  first  arrival  of  non-zero  values  of  a  solution  u(t,X)  of 

(6.1.2).   For  this  reason  we  refer  to  v    as  the  disturbance  speed.   Thus 

max 

(lU)  implies  that  the  component  of  the  group  velocity  in  the  direction  of 
the  normal  to  the  wave-front  does  not  exceed  the  disturbance  speed  of  the 
corresponding  non-dispersive  system.  [The  distiorbance  speeds  of  the  dis-      ■ 
persive  and  corresponding  non-dispersive  systems  are  probably  the  same. 
This  is  clearly  true  for  the  dispersive  differential  equations  (2.3-5)  because 
characteristic  hypersurfaces  and  normal  speeds  depend  only  on  the  "principal 
part",  i.e.  the  first -derivative  terms,  of  these  equations.] 
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7-   Electromagnetic  Waves  in  Isotropic  Media 
In  chapter  2  we  motivated  the  introduction  of  the  general  system 
of  equations  by  discussing  the  electromagnetic  field  equations  for  dis- 
persive media.   We  return  now  to  those  equations  in  order  to  see  how  our 
general  results  apply  to  them.   We  consider  here  only  the  case  of  isotropic 
media.   This  means  that  the  dielectric  and  magnetic  permeability  e'  and 'e' 
are  scalar  functions  of  cu  and  X.   Then  from  (2.1.24)  and  condition  1  we 
see  that  e"  =  £  -  (±X)~o     +   0(X~  )  and 


^^ 


^1^3 

0 

,»- 

'ih 

0 

^3  = 

0 

B2I3J 

0 

5^13] 

1 

0 

0" 

0 

1 

0 

0 

0 

1 

We  also  impose  conditions2-6  which  are  listed  in  appendix  G.   They  imply 
obvious  conditions  for  'e'^  and  e^.   Specific  examples  satisfying  these  con- 
distions  are  discussed  in  appendix  D. 


7.1  The  dispersion  equation  and  the  null  eigenvectors 

The  matrix  G  =  k  A  -  u^   may  be  obtained  from  the  above  ex- 


pression f or  (^  and  from  (2.1.21).   Thus 


G  = 


c(K) 


-c(K) 


-"^2^3 


(1) 


Let  us  denote  the  null  eigenvectors   of  G  by  r^  =   (R^,   R   )  where  R^  and  ^g 
are   3-vectors.      Then  the  equation  Gr_  =  0  becomes 

(2) 


W£-,R,    +  cK  X  R^   =  0, 


cK  X  R^    -  cje  R^   =  0. 
^      /^l  2""2 
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We  first  note  that  (2)  is  satisfied  with  w  =  0  and  KXR,  =KXR=0, 

i.e.,  oj  =  0  is  a  root  of  the  dispersion  relation  det  G  =  0.   However  for 

this  root  the  group  velicity  G^ =  (g  ,  g   ,   g  )  is  zero  because  g^  =  5w/Sk^  =  0; 

hence  the  rays  never  leave  the  "source  region"  (i.e. ,  the  projection 

of  the  manifold  v^  into  X-space).   Thus  the  zero  root  does  not  contribute 

to  the  asymptotic  solution  outside  of  the  source  region  and  is  therefore 

uninteresting. 

If  e^  (o))  vanishes  for  some  real  value  w  =  w  then  (2)  is  satisfied 
1^  o      ^ 

with  u)  =  0)  ,  R^  =  0,  and  R,  parallel  to  K.   For  such  solutions  the  group 

O'  ^2  —1  ^  ^  D     jr 

velocity  is  also  zero   (because  oj  is  independent  of  K) ,  and  the  corresponding 

asymptotic  solution  is  again  uninteresting.   We  shall  not  consider  either 

of  the  solutions  w  =  0  or  w  =  oj  further. 

o 

If  ojE-L^O,  (2)  implies  that 

R  -K  =  0  f3) 

and 

Ro  =  -^  K  X  R_  (k) 

The  last  equation  detemiines  R  in  terms  of  R  and  from  (3)  we  see  that 
there  are  two  linearly  independent  eigenvectors  r  =  (R  ,  R  ) ..  r  =  (R  .  R  ): 

ANA  ^/s^3_  ^       A^2       A^*  ^  "-^  '       A^2 

i.e.,  q  =  2.  If  we  insert  {\)  in  (2),  expand  the  vector  triple  product, 
and  use  (3)  we  obtain  (oj  e^e^  -  c  k  )R  =  0.  This  mplies  that  the  dis- 
persion relation  is  given  by 
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^^E^i^,X)e^{Lo,X)   =  c\^. 


The  matrix  A  is  given  by 


(jj 


^""^l^h 


("^2)03^3 


hence  the  orthonormality  condition  (3.2.U)  becomes 
(coe-,)  R^B:J  +  (oje^)  R^-R^  =  S.  .  • 

But  (4)  Implies  that 

j_   .    2  2    ..£-,.. 
(we^)  Eg 


Hence 


and 


where 


12    12 

R,  .R,  =  r:;.r^  =  0  , 

^2.  —1   —2  -^2  ' 


,j^2 


j^2 


£^(R^)"  =  e^CR^)"  =  ^ 


^^^^l^c    ^^^2)c 


"^1^2 


"1        "2 
We  also  note  that 

Thus  condition  7  is  satisfied  with 


{^E^)(^E^)^H^e^){^e^)^ 


^1   ^2 


(5) 


(6) 


(7) 


(8) 


(9) 
(10) 

(11) 


(12) 


(13) 
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From  the  dispersion  relation  (5)  it  is  clear  that  w  is  inde- 
pendent of  A  =  K/k.  The  function  oj  =  h(k,X)  is  defined  implicitly  by 
(5)  and  is  in  general  multiple -valued.  The  inverse  function  is  however 
single -valued  and  is  given  explicitly  by 


to 


k  =  m(u;X)  =  ^-  \/e^E^     .  (l^) 

It  is  defined  for  all  real  co  such  that  e  (cjj)ep(w)  >  0.   From  (ll)  and  (lU) 


it  is  easy  to  show  that 


^  =  TT^  (m'  =  m  ).  (15) 

^   2tan'  ^      w'  ^   ' 


Since  g^  =  -^  =  h'(k)a^,   a^  =  k^k,  we  see  that 

V 

G  =  (g^,  gg,  g^)  =  h' (k)A  =  A/m' (w) ,  (16) 


and 


g2  =  G.G  =  -^  =  ^  =  ^  .  (17) 
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7.2  The  polarization  vector 

Since  the  matrices  A  are  real  and  constant,  and  ^  is  real, 
ve  see  that  conditions  7?  8,  and  9  si's  satisfied  with  3  =  0  and  tj  given 
by  (7.1.13)-   Furthermore,  as  we  have  seen,  q  =  2.   Therefore  the  formulas 
(3.10.3^^^6)  for  z(t)  are  applicable.   However,  rather  than  use  these 
formulas,  we  prefer  to  return  to  section  3-9  and  derive  equivalent  formulas 
for  z,  in  which  the  physical  interpretation  is  emphasized. 

The  unit  vector  A  =  K/k  is  a  function  of  (t,X).   Let  A,  N^,  and 
B  be  mutually  orthogonal  unit  vector  functions  of  (t,X)  which  satisfy* 

B  =  A  X  N,  W  =  B  X  A,   A  =  N  X  B .  (l) 

1      2 
From  (7.1.3)  we  see  that  we  may  choose  any  two  directions  f or _R  and  R 

which  are  orthogonal  to  K  =  kA.   Therefore  we  set 


,1_,  Zip       p2  _    I 


Rt  =\/^B,      Rf^  =>/  ^  N.  (2) 

^1   V  e-j^'~'     --1   V  ^1  " 


Then  it  is  easily  seen  that 


and  (7.1. 3, U, 7)  are  satisfied.  From  (3.5.I)  we  note  that 


Therefore 


Z^  =  \/^  (ct  R  +  a  W),    Z^  =  \/^  (-OnW  +  0^B), 
'^l   V  £-,   1'"    2~  '    **2   V  Ep^   I"*"    2«' 


(5) 


*From  (7.1.16)  we  see  that  the  unit  vector  A^  is  tangent  to  (the  space 
projection  of)  a  ray.   Later  we  shall  take  W  and  B  to  be  unit  normal  and 
binormal  vectors  to  the  ray. 
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We  now  introduce  the  (complex)  polarization  vector, 

P  =  (p,B  +  p^N)  =  P'  +  1  P",  (6) 

where  P'  and  P"  are  vectors  with  real  components  and  B-,  ,  B_  are 
*"      "-  1  d 

defined  by  (3.9.6).   Then  (3.9. 12)  implies  that 

(p.)2  +(p..)2  ^p^.p  ^  |p^|2  ^  |p^|2  ^^^  (^) 

and  from  (l)  and  (5)  we  see  that 

Z,  =  /^  P,     Z^  =  /^  A  X  P.  (8) 

By  taking  the  real  and  imaginary  parts  of  (6)  we  find  that 

P'  =  3i  B  +  3^  N,     P"  =  p;[  B  +  p^  W  (9) 


where 


3^  =  p^  +  i  p|[,       P2  =  P^  +  i  Pg-  (10) 

It  is  clear  that  p'   p'  and  p"   p"  satisfy  (3. 9. 10),  therefore 


|P'|  =  \/(Pi)^  +  (Pp^  =PS    irl  =  /(PP^  +  (PP^  =P"'    (11) 
where  p'  and  p"  are  independent  of  t.   Thus 

P'  =  p'  Fn  cos  a'  +  B  sin  a'l  ,    P"  =  p"  [n  cos  a"  +  B  sin  a"  I  ,    (12) 

where  p'  sin  a'  =  P^ j   P'  cos  a'  =  3',  etc.,  and  from  (3. 9. 10)  (with  t 

replaced  by  x)  we  find  that 

t 
a'(T)  -a'(t)  =  a"(T)-a"(t)  =5  =  jT  T2(t')df.  (I3) 
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(7.2) 


The  real  and  imaginary  parts  P'  and  P"  of  the  polarization  vector  are 
illustrated  in  the  following  figure. 


The  lenths  p',  p"  of  the  vectors  and  the  angle  a' (t)  -  a"(t)  between  them 
are  fixed,  but  as  t  increases  the  pair  of  vectors  rotates  around  the  ray 
direction  A,  according  to  (13)-   [We  shall  shortly  investigate  the  geo- 
metrical significance  of  the  quantity  Xp  that  appears  in  (13).| 

From  (3.10.1)  we  recall  that  the  energy  density  w  which  appears 
in  (8)  is  given  by 


w(t)  =  w(t)  ^  ^^  \  -2/  n(t'  )<i^'  >  ' 

and  if  we  let  t  -> r  in  (3.7-5),  (^),  and  (8)  we  find  that 
w(t)  =  (z,A  z)  =  [(ce^)^  ZJ.Z^  +  (a3e2)^Z*.zJ, 

^   I — '      "^  "1 


and 


P'(t)  +  iP"(T)  =  P(t)  =  V^  Z,- 


(ih) 


(15) 
(16) 

(17) 


(7.2) 


Our  results  so  far  enable  us  to  determine  z(t)  =  (Z, ,  Z^)  along 
a  ray  in  tenns  of  the  value  of  z(t).   (Tlie  determination  of  z  is  discussed 
in  section  3-11  and  is  illustrated  in  chapter  5.)-  Let  us  summarize  these 
results:  From  (15)  we  may  determine  w.   Then  (l?)  and  (12)   (at  t  =  t) 
may  be  used  to  determine  the  constants  p' ,  p"  and  a' (t) ,  a"(T).   Once  t^ 
is  computed  (13)  yields  the  values  of  a' (t)  and  a"(t),  (12)  yields  the 
values  of  P(t)  =  P' (t)  +  iP"(t),  and  (8)  and  {ik)   provide  formulas  for 
■i  Z,(t)  and  Z-(t). 

It  remains  to  compute  the  value  of  t^  which,  from  (13),  determines 
the  rate  of  rotation  of  the  real  and  imaginary  parts  of  the  polarization 
vector  P:  From  (7.1.6,9)  we  find  that 

(r^,A°  r^)  =  0  (18) 

therefore  (3.9.^)  yields 

Tp  =  iL>   A°r2)  +!_(/.  A%^  ).  (15) 

v=l  V 

But  from  (2.1.21),  (2),  (3),  and  (7-1.17) 


y~  (r^,  aV  )  =  crR^.V  X  R?  -  r|.V  x  R^I  =  -  / --^   [w-V  x  W  +  B-V  x  b] 
v=l        ^v       '^  '^      V   12 


1 
2 


g  fw-V  X  W  +  B.V  X  b1.  (20) 


By  using  (E.I3)  of  appendix  E,  we  now  find  that 


(r-^,  A^  r^  )  =  -g  N.(A.V)B  =  -N- (G.V)B.  (21) 

V=l  V 

From  (7.1.6),  {2),    (3)  we  see  that  the  first  term  of  (19)  is 

^i'   A°r2)  =  (-e^)4-(R^)^+(-e2)^Rj.(R2)^  =  o^B-N^  -  a^N'g^,       (22) 


(7.2) 
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where  <x  =  {'^/e   ){'^s   )  .   From  (7.1.11)  we  see  that  ou  +  ap  =  1,  therefore 

(r-^,A°rf)  =  ch(B.W,  +  K-B,  )  -  N-B,  =  cu  (B-W) , -N-B,  =  -N-B,.         (23) 
By  inserting  (23)  and  (21)  in  (I9)  we  obtain 


5t*«-' 


dB 

-I'di  •  (211) 


In  {2k)    d/dt  denotes  differentiation,  with  respect  to  time, 

along  a  ray.   We  now  choose  N  and  B  to  be  unit  normal  and  binormal  vectors 

to  the  projection  of  that  ray,  which  is  a  space -cur ve .   Since  A  is  a  unit 

tangent  vector  to  the  projection.  A,  N,  and  B  satisfy  the  Frenet  equations 

dA  dU  dB 

dS  =  ^  H'    d?  =  -'^  A  +  T^B,   ^  =  -T^N.  (25) 

Here  k   is  the  ciurvature  and  1     is  the  torsion  of  the  projected  ray;  o   denotes 
an  arclength  parameter.   But  (7.I.I6)  and  the  ray  equations  imply  that 

-£  =  G  =  ^k  =  ^^,  (26) 

dt   -^   m'  '-   m'  da  '  ''   ' 

hence  dt  =  m'  dcr.   Thus  {2k)    and  (25)  yield 

-j_     dB   T 
-^2  =  -  J^  S>-  ^  =  ^  ^  (27) 

and  from  (13) ^ 

t  a(t) 

5  =  J  ^  df  =  J        x^  da.  (28) 

T  a{'x) 

Thus  we  see  that  the  rate  of  rotation  of  the  vectors  P'  and  P"  is  de- 

termined  by  the  torsion  t  of  the  projected  ray.   If  this  space  curve 

is  a  straight  line  (as  it  is  for  a  homogeneous  mediiim)  t  is  zero  and  the 

vectors  P'  and  P"  are  constant  on  the  ray.   If  the  projected  curve  is 

a  plane  curve,  x  is  again  zero  and  the  vectors  P'  and  P"  remain  fixed 

in  the  A,  N,  B-frame. 


(A) 
Appendices 
A.   A  theorem  for  homogeneous  media 

Theorem.   If  A  (v=l,...,n)  and (^  are  independent  of  X  and 
condition  6  holds  then 

T   =  ri  V.G  +  T1.1&.  . 
to   '-2   '^   'I-'    to 

n 

Proof ;        Our  assumptions  imply  tliat  G  =/_  k  A   -a)0('^)      is 

v=l 

independent  of  X.   Therefore  the  roots  oa  =  h(K)  of  det  G  =  0  are  independent 
of  X,  and  it  follows  from  the  ray  equations  (3'3-l)  that  k^,  w,   and  G  are 
constant  on  each  ray.   Thus  the  null  eigenvectors  r  of  G  may  be  so  chosen 

that  they  are  constant  on  each  ray,  and 

n 

Zm    d   m   _  -,  /  -,  N 

gr   =  —  r   =  0;  m=l,...,q.  (1) 

v=0      V 


^.    ^     .o,      Sco    Sh     _     „,l«o     „ 

Since  G  =  -A  and  g.=^r  =^rj   G,   =A'^-Ag..   Hence 
u)          °j   dk .   dk .     k .  J 

J     J      J 

^ ,m  =  1 , . . .q: 
(r^  A^r"^  )  =  (rS  G,  r"^  )  +  (rS  g.A°r"^  );     .   ,  (2) 

J  J   J  J 

From  (3.5.6  ,7)  and  (2) 

n 

1/  i  ,0  m^    _     /  t  .0  m  V  ,  \  f    I      . j  m   s 
T.   -  ;r(r  ,A,r  )-  ti.5.   =  (r  ,A  r   )  +  /_  (r  ,  A""  r   ) 

o    J=l  J 

But  from  (l)  we  see  that  the  last  term  in  (3)  vanishes.   Furthermore,  since 

G,  r^  +  Gr^  =  (Gr-^),   =  0, 
k.-^     -k.   ^  '~  'k.    ' 
J        0        J 


Hence  (3)  and  {k)   yield 
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(A) 


Since  k.  =  s   is  constant  on  a  ray 


Sk.   dk. 


n  n 

J  O     V=l       J  V     V=0         V 


0, 


(6) 


therefore 


V=l    V    ""   j  v=l 


n 


=  7_  A  s    .   (7) 
X  X.   ^—x  x,x. 

V  J   v=0      V  J 


Furthermore 


n 


\    m   j_   \ 


r   =  /    r,    x~^  =  /   s     1^1  - 
-x    4-r-     -k.   CK,^   ^  XX  -k 


Wow  (3-5.6),  (8)  and  (?)  yield 
1/  t   „o  m 


n 


-^.-|(£.^*t£")-\^.  =  ^   (i^*'*"^;' 


(8) 


n    n 


-   I_  >_  irj,   A^  s 
v=0  j=l 


m  - 
x.x,,  i^k.' 


n 


t  „   ..m  \     S       '  ^       -  "1 


\ 


(9) 


=  -/_  (G  r"  ,  r"'  ) 

(The  last  equation  in  (9)  follows  from  the  fact  that  G  is  hermitian. ) 
We  now  compare  (5)  and  (9)  and  find  that 

"^to  "  '^mt  :  (10) 

Then  (10),  (3-7-3)j  condition  6  and  (3. 2. 11)  yield  the  required  result, 


T.  =\\y,      +t''J  =  I  (V.G)S,   +  (rS  c^^r"")  =  [i  V-G  +n  l6„  .   (ll) 


'im   2  ■-  to    mt-'    2 


to 


'-2   -    't-J  em 
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B.   Generalized  energy  density  and  Poynting  vector 

If  u  Is  any  real  solution  of  the  general  system  of  equations 
(2.1.26,  7)  we  define  the  generalized  Poynting  vector  S  =  (S,  ,  ...,  S  ) 
corresponding  to  u  by 

S^  =  2(u,  A^  u);     V  =  l,...,n.  (l) 

For  real  solutions  u  =  (E,  H)  of  the  electromagnetic  field  equations  we 
see  from  (2.1.21)  that  if  K  is  an  arbitrary  vector, 

K-S  =  >   k  S  =  2crH.K  X  E  -  E-K  X  h]  =  llcK-E  X  H       (2)    . 
v=l 

hence 

S  =  l6n  S^,   where  S^,  =  ^  E  X  H.  (3) 

But  &^   is  the  usual  "Poynting  vector"  of  electromagnetic  theory.   For 
this  reason  we  have  called  S^ the  "generalized  Poynting  vector". 

For  complex  solutions  we  define  S  in  terms  of  the  real  part  of 
u.   Thus 

S^  =  2(Re  u,  A^  Re  u)  =  |[(u  +  u*)  ,  A^  (u  +  u*)J 

=  |[(u,a''u)+(u*,a',u*)+(u*,a"u)+(u,aV)]  .  (k) 

Throughout  this  report  we  have  studied  complex  asymptotic  solutions  of 

iXs 
the  form  e    z.   For  such  solutions  (k)   becomes 

o    If/   A^  ^  /  -x-  ,^  •x-\  .   2iXs/  j(.  .V  ,     -2iXs/   .v  ^,-1      ,  , 
%   =  i[(2,A  z)+(z*,A  z  )  +  e    {z*,   A  z)  +  e     (z,A  z  )J  .     (5) 

We  note  that  the  last  two  terms  in  (5)  contain  oscillatory  factors,  hence 
if  we  average  (5)  over  a  short  time  interval-^,  the  average  of  the  last  two 


*A  space -average  or  space -time  average  will  yield  the  same  result. 


(B) 
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terms  will  be  asymptotically  zero  for  X  ->».   Thus  the  average  value  of 
S  is  given  by 

<S^>=|[(zyz)  -f(z-,  aV)].  (6) 

If  we  assume  that  A  is  real  for  v  =  l,...,n;  then 

(z*,  aV)  =  (z^A^'z)*  =  (A^'z^z)  =  (z^a'^z)  and  (6)  becomes 

<S^>   =  (z^z).  (7) 

The  vector  <S  >  =  «S  >,...,  <S  >)  is  the  average  Foynting  vector  of  the 
given  asymptotic  solution. 

¥e  now  introduce  (3. 5.1)  in  .(7)  and  we  use  the  basic  identity 
(3.2.8)  and  the  definition  (3-7.5)  of  w.   We  find  that 

hence 

<S>  =  wG.  (9) 

We  also  find,  from  (3-7.7)  that 

Thus  if  our  system  is  conservative  and  the  matrices  A  ,  ...,  A  are  constant, 
then  Tj  =  3  =  0  and 

w^  +  V.<S>  =  0.  (11) 

This  equation  is  the  well  known  energy-conservation  equation.   It  explains 
why  we  have  called  w  the  "energy  density".   To  be  more  precise  we  should 
call  it  the  average  energy  density.   From  (9)  we  note  that  the  Foynting 
vector  <s)  has  the  direction  of  the  (projected)  ray  and  the  magnitude  wg, 
where  g  is  the  group  speed.   Thus  it  is  particularly  clear  in  the  asymptotic 
theory  that  the  Foynting  vector  measures  the  energy  flux.   From  (3)  we 


(B) 
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see  that  the  average  energy  density  of  the  electromagnetic  field  is  given 
by 

w^  =  w/l6n.  (12) 

It  is  interesting  to  note  that  equations  1,  h,    5,  6,   7 ,   8  remain 

valid  with  v  =  0,  hence  v  =  (S  )  where  S  =  2(u,  A  u) .   Now  the  matrix 

A  (w)  and  the  above  quantity  S  are  defined  only  for  asymptotic  solutions, 

because  then  w  =  -s^  is  defined.   For  exact  solutions  S  is  -undefined 

t  o 


CJ. 


unless  the  system  is  non -dispersive,  for  then  A  = <5   is  independent  of 
For  the  non-dispersive  electromagnetic  field,  it  is  easy  to  show  that 
S  =  l6ic  S  „,  where  S  _  is  the  usual  electromagnetic  energy  density, 

O  Oilj  Oilj 

^oE  =  fe  t^'^^^'j-  ^^3) 

For  dispersive  systems  it  does  not  appear  to  be  possible  to  give  a 
reasonable  definition  of  the  energy  density  in  terms  of  exact  solutions, 
but  only  in  the  framework  of  the  asymptotic  theory.   A  similar  remark 
applies  to  the  group  velocity  vector,  G_. 

As  we  have  seen  in  section  (3. 11)  an  asymptotic  solution  is  often 
given  by  a  sum  of  terms  of  the  form 

iXS/ 


(j) 


u(t,X)  =  '^  u,  ..;  ii/  -N  =  e       z  /.>, 


(1^) 


(Each  ray  that  passes  through  the  point  (t,X)  contributes  one  temi  to 
this  sum.)  For  such  solutions,  {h)   becomes 


(B) 
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S     =|r     E\e      ^     (j)      ^"^)^    {z..,,a\,    J-te      ^^J)      («^)^    (z^,,aV   J 

iX|s/.^+S/    n  ixf-S/.v-S/    vl  ^ 

+  e      ^(J)      (-)      (z-..,a\,    Jfe      ^     (j)      (-)      (z^.yz!   J)  (15) 

We  assume  that    |s.|    ^    |s    |    for   j   ^  m.      Then  by  averaging   (I5)   we   find 
that 

Thus 

<S>  =Z  <S(.)>  (IT) 

J 

where  (S/.\)  is  the  average  Poynting  vector  for  the  asymptotic  solution 

iXs,  .V 
U/  .N  =  e   ^   ^^-V  '^■^  ^f  •)    ^^   *^^  energy  density  of  the  same  solution, 

then  (11)  Implies  that 

Thus  we  see  from  (l?)  and  (I8)  that 


^t  +  ^'^^^  =  °  ^^5) 


where 


w  =  X!!  w^-v  (20) 

From  (19)  and  (20)  we  conclude  that  the  energy  density  of  an  asymptotic 
solution  consisting  of  a  sum  of  terms  (l4)  is  the  sum  of  the  energy  densities 
associated  with  each  term. 
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(c) 


C.   Solution  of  the  equation  for  the  energy  density 

The  energy  density  function  w  satisfies  the  equation  (3.7-7) 
^   +  [V.G+  2(Re  ^+   3)]w  =  0.  (l) 

In  order  to  solve  this  equation,  we  set 

r    r  \ 

V  =  Jar,  7  =  exp  (^  -2  /[ReT(t' )+p(t' )]df  >  ;  (2) 

^^  *o  ^ 

then  it  is  easily  seen  that  a  satisfies  the  simpler  equation, 

i  -^  (^-2)^  =  0.  (3) 

We  nov  set  t  =  x  ,  s  =1  and  introduce  the  (n+l) -vectors 
o'   o 


§  =  (g^,  g3_,  -..,  g^)  =  (g^,G) 


^  -  '  ax  '    ^yi'    '•'    '    ^x    j  "  i  ax  '   ^ 

O       1  V        ^    O 


(h) 


^   .   ^.        d     -^  .^  ^.         da 


Then  from  (3) 

V-(a2)  =  Va-G  +  a(V.G)  =  ^  +  a(V-G)  =  0.  (5) 

Let  »  be  an  arbitrary  point  set  in  the  space  (jot   parameters 
r  which  label  the  rays  (t,X)  =  [*  jS.^''^»r)]  •   Let  J  be  the  set  of  points 
in  space -time  filled  by  the  "tube"  of  ray  segments  for  which  r  is  in  ^ 
and  t  <  t  <  t  .   The  region  O  is  illustrated  in  the  following  figure: 


't 
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We  now  apply  Gauss'  theorem  to  the  region  O 
follows  that 


From  (5)  it 


0  =  j  V.  (oG)  dx^dx^  ...  dx^  =  j     oG-W  dS, 


(6) 


rays  G-W  =  0.   On  the  hyperplanes  t  =  t  and 


where  W  is  the  unit  outward  normal  vector  to  the  surface  lA     of  o    .      On  that 

portion  of  s^   generated  by 

t  =  t   G.N  =  +  g  =  +  1.   Hence  (6)  becomes 

it^  t^j 

But  dx  . . .  dx  =  j(t;r)d7^ . . .d7  =  j(t)dr,  therefore 


■  a  dx-,  . .  .  dx  =0. 
n 


(7) 


^  [aCt^) jCtg)  -  a(t^)j(tJdX  =  0, 


(8) 


and  since  ^  ,   t  and  t  were  arbitrary  we  conclude  that  a(t)j(t)  is 

constant  on  a  ray,  or 

J(t  ) 


a(t)  =  a(t^) 


o'  j(t) 


(9) 
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(C) 


Wow  (2)  implies  that  ■w(t  )  =  a(t  ),  therefore  from  (2)  and  (9)  we  find 


that  the  solution  of  (l)  is  given  by 

J^*o)       f     r 

w(t)  =  w(tj  -r7%v  exp(  -2      [ 


(t)  =  w(t^)  ^^  exp/  -2j  [ReTl(t•)+3(t•)]dt• 


(lO) 


In  order  to  solve  equation  (3.6.2)  we  set  w  =  0  and  p  =  0. 
From  (1)  we  find  that 

if  +  [|^'5  +  ^]°  =  0.  (11) 

while  (10)  yields 


a(t)  =  a(t^) 


J(t^) 


j(t) 


.1/2        r    t 

exp  <  -/   Ti(t')dt' 


(12) 


Since  (3-6.2)  and  (ll)  are  the  same,  (12)  is  the  required  solution. 
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D.   The  dielectric  permeability  function 

Theoretical  derivations  LlU  J  of  the  dielectric  permeability 
for  isotropic  media  with  a  single  "resonance  frequency"  lead  to  formulas 
of  the  form 

Bi(-.X)  -  1  =  ,,f  ,  .  (1) 

p  -co  -ivoj 

Here,  as  in  chapter  2  and  chapter  7j  ^  denotes  the  frequency  and  £-. 
denotes  the  dielectric  permeability,   (j)  is  proportional  to  the  number 
density  of  oscillators  of  resonance  frequency  p,  and  v  is  proportional  to 
the  damping  force.  ^,  p,   and  v  may  be  furnctions  of  X,  and  v  >  0.   For 
media  with  n  resonance  frequencies,  the  right  side  of  (l)  is  replaced  by 
a  sum  of  n  such  terms.   For  an  isotropic  plasma,  (l)  is  valid  with  p  =  0. 
Then  (j)  is  called  the  "plasma  frequency"  and  v  the  "collision  frequency". 

If  (f)  is  independent  of  X,  we  set  X  =  (|).   Otherwise  we  take  X 
to  be  some  average  value  of  (f)  (or  an  average  value  of  p)  and  we  set 
r  =  p/x,  p  =  (f)/x.   Then 

p  =  Xr(X),   (f)  =  Xp(X).  (2) 

As  in  section  2.1  we  set 

oj  =  w/X;    £^(^,X)  =  e\(J^,X).  (3) 

Then 

r  -co  -lOJv/X 
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^- 


where 


¥e  now  expand  {h)    for  large  X  and  find  that 
e^  =  £^  -  (iX)"^6^  +  0(X"^) 


^1  =  1  - 


P  ., 

2  2  ' 


wvp 


CO  -r 


1  ,   2  2.2 
(w  -r  ) 


(5) 


(6) 
If  we  take  the  magnetic  permeability  of  the  medium  to  be  e     =   Ij,  then  from 


(2.1.21+) 


s,l3    0 


3^ 


=  <J-(ix)"^<=^+  o(x"^),   (f= 


£,13   0 


^  = 


6,13   0 


(7) 


Thus  we  see  that  condition  1  of  appendix  G  is  satisfied.   From  (6)  we  find  that 

2,  2   2^ 

(8) 


a^  ("  ^1^  =  1  ^   /2   2.2   ^ 


(co  -r 


■N  ' ' 

00^  ^ 

hence  A  =  '^  {^C>)   >   I/-^  where  I^  is  the  6X6  identity  matrix.   Since  lira  G  =  I^ 

w  — >  00 

we  see  that  condition  11  is  satisfied  with  A  =  I^.   Condition  2  then  follows 
from  condition  11.   Condition  k   is  easily  verified  frora  (7)  and  (k) ,   and 
condition  3  may  be  verified  by  direct  computation. 

In  order  to  verify  condition  5  we  note  that  the  equation  det  G  =  0 
is  given  by  (7.1. 5)  with  e^  replaced  by  e,  and  e  replaced  by  1.   Thus  (k) 

yields 

2  2    ,  2,  2    2v  ,  2    2    .   /,  >  ,^. 

cjp  =(ck   -oj)(r   -co   -  icov/X).  (9) 

It  is  necessary  to  show  that  the  roots  of  (9)  lie  in  the  half -plane  Im  w  <  0. 

For  cj  =  00  it  is  easy  to  find  the  roots  and  they  are  all  real.   For  large  X 

a  perturbation  solution  of  (9)  shows  that  Im  co  <  0.   Since  the  roots  of  (9) 

cannot  be  real  for  any  finite  X,  the  roots  must  remain  in  the  lower  half- 

plane  as  X  decreases.   The  remaining  condition  of  appendix  G  are  all  easily 

verified. 
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E.   Some  vector  identities 

¥e  begin  with  the  standard  vector  identity 

Vx  (A^XAJ  =  A^(V.A_)-A^(V.A^)+(A  .V)A_-(A^.V)A_.  (l) 

If  A,,  A^,  A_  are  orthonormal  then  (l)  yields 

A»»-L    l^d  •~j 

A,  .V  X  (A_XAJ  =  A^  •  [(A  -VU-CA^-VUJ  .  (2) 

We  now  apply  (2)  to  the  vectors  A.  W,  B  which  satisfy  (7.2.1).  This  yields 
N.V  X  N  =  K-V  X  (B  X  A)  ="  N-  r(A-V)B  -(B.V)aI  ,  (3) 

B-V  X  B  =  B-V  X  (A  X  W)  =  B- [(N.V)A  -(A-V)n] ,  [h] 

A-V  X  A  =  A-V  X  (N  X  B)  =  A-  r(B-V)W  -(W-V)bI  .  (5) 


and 
But 


W-(A-V)B  +  B-(A.V)W  =  (A.V)(B-N)  =  0,  (6) 

A- (B-V)N  +  W.(B-V)A  =  (B-V)(N-A)  =  0,  (?) 


«w*   •^*'    '^v^ 


and 


B-(N-V)A  +A-(N-V)B  =  (N-V)(B-A)  =  0.  (8) 


If  we  add  (3)  and  {W)    and  subtract  (5)  we  find  that 

W-V  X  W  +  B-V  X  B  -  A-V  X  A  =  2N-  (A-V)B.  (9) 

Wow,  by  definition,   K  =  kA,  hence 

VxK  =  VkXA+kVxA  (lO) 

^  -~  -^ 

and 

A-Vx  K  =  kA.V  X  A.  (11) 

Therefore,    since  K  =  Vs, 

A-V  X  A  =   (l/k)A-V  X  K  =  r  A-V  x  Vs  =  0,  (12) 

and   (9)   becomes 

W-V  x  W  +  B-V  X  B  =  2N- (A-V)B.  (13) 


-100- 

(F) 
F.   Notation 

In  the  system  of  equations  studied  in  this  paper, the  independent 
variables  are  t  and  X.   X  is  a  vector  with  n  components.   Tlie  dependent 
variables  u  and  v  are  column  vectors  with  m  components.   In  general, 

A  denotes  an  n-vector  with  components  A,,..., A  or  a  ,...,a  .   If  n  =  3, 
vw  _L      n      J-      n 

A  is  an  ordinary  3-"vector. 

a  denotes  a  column  vector  with  m  components.   To  save  space  we  often  write 

the  components  horizontally  rather  then  vertically.   Thus 

a  =  (a.,  , .  .  .  ,a  ) . 
'^     1      m 

A  denotes  an  m  x  m  matrix  or  a  scalar.   Thus  Aa  is  a  colimin  vector, 
a  denotes  a  scalar  or  (rarely)  an  n  x  n  matrix. 

The  inner  product  of  two  column  vectors  a  and  b^  is  a  scalar 
defined  by 

j=l  ^   ^ 
The  star  denotes  the  complex  conjugate  (or  the  Hermitian  conjugate  for 
matrices).   The  following  properties  of  the  inner  product  are  used.   Tliey 
follow  easily  from  the  definition 

(a.b)*  =  (b.a);   (ca,b)  =  c^(a,b);    (a,cb)  =  c(a,b); 
(a,Ab)  =  (A*a,,b);   if  A  is  hermitian  (a,Ab)  =  (Aa.,b). 

The  summation  convention  with  respect  to  repeated  indices  is  used. 
Usually  the  indices  run  from  1  to  n.   If  not,  a  parenthetical  note  is 
inserted.   Thus 


a  b    denotes  x   a  b 
v=l 


denotes   /   a .b . 

Some  column  vectors  a  with  6  components  are  defined  by  an  ordered  pair 
of  two  3-vectors.   Thus  a  =  (A,B)  =  (A-^^A^^A^^B^jB^^B^) .   Tlie  symbol  (A,B) 
is  not  an  inner  product.   Only  column  vectors  a,  b  appear  in  inner  products. 
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G.   Summary  of  Conditions 

The  following  basic  conditions  are  imposed  throughout  this  paper: 

1.  (|'(w,X;A)  =  (^(t^,X)  -  (iX)""^(>Q(oj,X)+0(x"^).   For  real  co,  G    is  hermitian. 

2.  For  real  co,  A  (lo,X)  =  -^  (loS)  is  positive -definite. 

3.  For  real  u  and  for  K  =  (k  , ...,  k  )  with  real  components,  the  multiplicity 
of  each  eigenvalue  of  the  matrix  G  =  k  A  (X) -ojfe  (lo,^}  is  independent  of  oj  and  K. 
^'     G  ('^jX;X)  is  an  analytic  function  of  w  in  the  half -plane  Im  w  >  0. 

5.  For  K  =  (k^,...,k  )  with  real  components  the  solutions  u)  =  h(K,X;X)  of  the 
determinantal  equation  det  G  =  det  |_k  A  (X^ -co  G  (lo,X;X)  |  =  0  lie  in  the  half- 
plane  Im  to  <  0. 

Conditions  1,2,  and  3  are  used  repeatedly  in  the  expansion  procedure  of  chapter 
3  and  the  perturbation  analysis  of  chapter  h.      Conditions  k   and  5  are  used  in 
obtaining  the  fourier  integral  representation  of  chapter  k   and  its  residue 
evaluation. 

The  following  special  conditions  are  required  for  certain  purposes.^ 
(r  , .  .  .  ,  r;  are  null  eigenvectors  of  G)  : 

6.  For  real  co,  oO  is  hermitian. 

7.  [r^,  u^(co)^J]  =n('^)5.  .;     i,j  =l,...,q. 

n 

8.  -  /   (r^,A^  r"^)  =  2p  6,  ;    t,m  =  l,...,q. 

v=l  V 

9.  The  hermitian  matrices  A  (X),...,A  (X)  and   g(w,X)  are  real  (hence 
symmetric)  for  real  u. 

10.  For  real  <^  ^  0,      §{-bi,X)    =  ,f  (w,X)  and  <7C^(-co,X)  =  -  25'(oj,X). 
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11,  g(wjX,;X)  tends  to  a  positive  definite  matrix  A(X)  as  w  tends  to  infinity 

in  the  upper  half -plane,  and  for  real  w,A  (lo,X)  >  A(X) ,  i.e.  the  matrix  A  -A 

is  non-negative. 

Conditions  6,7,8,  and  9  are  required  for  the  solution  of  the  transport  equations 

in  the  various  cases  discussed  in  sections  3- 6-3-9 • 

Condition  10  is  used  in  section  h.'^   to  simplify  the  solution  of  the  system 

of  equations  with  constant  coefficients.   Condition  11  is  used  in  section  3 

to  derive  the  inequality  (6. 3.1^)- 
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